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A whole new family of Speedomax® Recorders 
has come to the aid of time-pressed scientists in 


the past few years . . . some units only within 
the past months. Here are a few of the many 
jobs they can do: 

If you measure temperature by means of re- 
sistance thermometry, for instance, the correct 
Speedomax gives you a 5 place readout; not as 
accurate as the L&N Mueller bridge’s 7-place 
readout, but ample for a variety of applications. 

An L&N Wenner potentiometer gives you a 6 
place readout in thermocouple thermometry, but 
the precision of an automatic Speedomax with a 
4 place readout may easily fill your requirements. 

In the measurement of pH, our automatic re- 
corder is comparable in accuracy with our high 
quality manual instruments. 

Some of the newest laboratory “right hands” 
are engineered time-savers. The X-Y Speedomax 
can plot one quantity against another (such as 
temperature and temperature difference) and save 


could 


3 individual plottings. With the X:-X2 Speedo- 
max, any two measurements can be plotted simul- 
taneously against time. To record power-level 
measurements in a matter of minutes instead of 
hours or days, the automatic Speedomax power 
level recorder is the answer. 

All of these Speedomax Recorders, and many 
others, are standard instruments. Even in the 
rare cases where there’s no standard recorder to 
meet a need, modifications are usually available to 
enable you to record any electrical signal. 

We can give impartial service or, if you wish, 
advice, on either manual or automatic instruments, 
for we produce both. You can get the facts from 
our nearest office or from 4992 Stenton Ave., 
Philadelphia 44, Pa. 
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Which speaker is making the sound? In echoless chamber at Bell Labs, Robert Hanson 
measures test subject’s ability to localize sounds—observes how two ears operate in 
partnership. This and other tests may point the way to better telephone instruments. 


In listening to stereophonic music, how 
is it that our ears and brain construct 
a picture of the entire orchestra with 
but two samples (the sounds from two 
speakers ) to work with? 


How is it that our ears and brain are 
able to pinpoint one voice in a roomful 
of talkers—to listen to it alone and 
ignore the rest? 

What makes two ears better than one? 

Bell Telephone Laboratories scientists 
are searching for the answers. For, in 
finding them, better telephone instru- 
ments and better ways of transmitting 
sound will surely result. 

Our hearing performs feats that no 
electronic system can yet duplicate. 
How? Laboratories scientists believe the 
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secret lies in the way our two ears func- 
tion in partnership and in the way our 
neural network connects them with our 
brain. The problem: to discover what 
functions the network performs and to 
see whether electronic duplication might 
enhance understanding. 

The work is under way. Electronic 
circuits that simulate the operation of 
nerve cells have already been created 
—and conceptual models of the neural 
network are being constructed. 

Alexander Graham Bell’s interest in 
deafness and hearing led to the inven- 
tion of the telephone. Bell Laboratories’ 
current explorations in binaural sound 
may well lead to important new ad- 
vances in the transmission of speech and 
music. 
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THEORY OF THE A. F. JOFFE METHOD FOR RAPID MEASUREMENT 
OF THE THERMAL CONDUCTIVITY OF SOLIDS* 


BY 
G. SWANN! 


W. F. 


ABSTRACT 


The figure below represents a block B of high thermal conductivity fastened to a 
specimen S of length / parallel to the axis of x, and of lower conductivity k. Att = 
all is at constant temperature. Then, suddenly, the right hand side of S is brought 
to a temperature r.. and kept there. If @ is the excess of the temperature at some 
point over r., the general problem consists in finding @ as a function of x and ¢. 

The simple theory assumes that 6, at x = 0, decreases so slowly with time that 
the temperature gradient in S is, at each instant, independent of x and equal to 6/1. 
This assumption leads to @ = 0o¢~**‘, where @p is the initial value of 6, and A, = k/Msl, 
M and s being, respectively, the total mass per unit area of B projected on the plane 
of contact, and the specific heat of block B. 

The more exact theory gives 


6= A,e7 et [cos _ 
0 l 


where the constants A, are to be determined to give 6 = 0) at t = 0. X, is given by 
R\,/\a = v2 where R is the ratio of the thermal capacity of the specimen to that of 
B, fy, is a solution of the equation Gy, tan Y, = 1, and G, is a constant function of s 
whose form depends upon the nature of the boundary conditions at x = 0—the place of 
junction of Band S. For the case where there is perfect thermal contact at x = 0, 
G, = 1/R and is independent of s._ For the case where there is thermal resistance at 


x = 0,G, = [R + ky,2/hl]" where h is defined by 


hA@ = heat transferred across the boundary per cm? per sec. 


and A@ is the temperature drop. It turns out that, at the junction, 6 is of the form 
@ = Age! + Ase ™! + Aye! + 


where Aj, Az, Az, etc. are increasingly larger than Xo, and by very substantial amounts, 


= Since this paper was completed, there came to the author’ s attention an unpublished 
copy (a preprint) of a paper by M. A. Kaganov dealing with the theory of Joffe’s method. 

1 Director, Bartol Research Foundation of The Franklin Institute, Swarthmore, Pa., and 
a Senior Staff Advisor for The Franklin Institute Laboratories for Research and Development. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JouRNAL.) 
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so that, except near ¢ = 0, the term involving Ao is all important, and 6 = Age! at 
x = 0. It then results that, for R small, and h reasonably large 


She R for thermal junction resistance absent 
Aa 
1 
3 
and 
1 
=- R for thermal resistance included. 
1 
* 3 


Xo is the experimentally measured value, so that we should have, for the two cases, 


k R 
+ 3 )r 
k re 
Mst ~ (1 ~gt 
The elementary theory omits the quantities in parentheses. 


The actual determination of the A,'s to complete the expression for @ as a func- 
tion of x and ¢ is carried out for the case where thermal resistance at the junction is 


respectively, 


neglected. 


A 


INTRODUCTION 


A method used extensively for measuring the thermal conductivity 


of solids, including semi-conductors, in connection with thermoelectric 
research, is dependent upon the following principle 2 


2A. F. Jorre, Can. J. Phys., Vol. 34, p. 1342 (1956). 
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A block B of high thermal conductivity is fastened to a slab S repre- 
sented by the shaded portion of the figure. The thermal conductivity 
of S is low compared with that of the mass B, so that we may assume, 
for B, a temperature which, though time dependent, is constant through- 
out the mass although the temperature may vary throughout S in a 
manner determined by the nature of the problem. 


Our problem is such that, at ¢ = 0, everything is at the constant 
temperature ro. Then, suddenly, (at ¢ = 0) the face AB of Sis brought 
to the temperature r,, in actuality the temperature of melting ice, and 
kept there for all time. We desire to know the temperature 7 of S at 
each subsequent instant, as a function of «, measured from left to right, 
the plane of junction with B being at x = 0. We assume the dimensions 
perpendicular to the axis of x to be so large that we may treat the prob- 
lem as one-dimensional. Moreover, we assume that B loses heat only 
at the plane of junction with S. 
We make the following definitions: 


M = mass of B divided by the area of its junction with S, 
s = specific heat of the material of the mass B, 
/ = thickness of slab S, 
k = thermal conductivity of the material of the slab S, 
o = specific heat of the material of the slab S, 
a? = ratio of thermal conductivity to product of specific heat and 
density for S, 
A@ = drop in temperature at the junction of S and B on account of a 
finite thermal resistance at that junction, and 
h = coefficient of A@ in the expression for the heat conducted per cm? 
per second across the junction. 


The solution of this problem obviously provides a method for de- 
termining the thermal conductivity of S in terms of the temperature- 
time relationship for B, and of the relevant dimensions and other 
parameters of the system. 

The basic differential equation controlling the temperature 7 in the 
slab S is 


Or 


= a? (1) 


It will be convenient to define a quantity @ by 


Te 


so that 6 is the excess of r over 7,. In terms of 6, Eq. 1 becomes 


at 


= a? V%, (2) 


i 
|_| 
ot ; 
q 
|_| 
2 
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Conditions Governing the Solution 
We have: 


A. @ must satisfy (2) for all values of t over the range 0 < x <1 
B. 6 = 6) at t = O for all values of x between x = 0 and x = 1 
C. At x = 1, 6 = 0 for all time. 

D. Att = «, 6 = 0 for all values of x, from x = Otox =1 

E. Atx =0 


(3) 


00 =Msk 0°60 00 
Ms — —— = k— 3.1 
according as we do not or do take account of a thermal resistance be- 
tween the junction of S and B. 


Customarily Used Approximate Solution 


If we assume that Ms is very large, it is clear that, unless k is very 
large, 6, at x = 0, will diminish very slowly with the time. If we as- 
sume that, at each instant, 00/dx assumes, in S, a value independent of 
x and equal to —6/1, Eq. 3, neglecting thermal resistance at the junction, 
yields‘ 

06 ké 


Ms 


= e—kt/ Mat (4) 
which expression enables k/ Ms to be determined through knowledge of 


the experimentally measured 6@,t relationship. 


3 In terms of the notation in the system of definitions, we have 
20 + 40) _ , 
at Ox 


where @ is the temperature at x = 0, so that the temperature of the surface of B at x = 0 is 
6+ A6. From these relations we deduce immediately 


06 0% 00 
Ms— — —— 
at h Ox 

‘It is clear that 00/dx must, in actuality, always differ from the constant value — 6/1 to 
the extent necessary to provide for the fact that however slowly 6 diminishes with time at 
x = 0, heat must flow out of the various parts of S at a corresponding rate sufficient to insure 
the continued fall in the average temperature of S. All of the heat flowing out of S on this 
account passes through x = J, while none of it passes through at x = 0. It is manifest there- 
fore that d@/dx at x = / must differ from that at x = 0 to the extent necessary to harmonize 
matters. 


SWANN [J. F. 1. 
00 a0 

Ms— =k | | 

or 
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THE GENERAL THEORY 


In view of the uncertainty as to the magnitudes of the errors in- 
volved in the assumptions leading to (4), and of the general applicability 
of the method, it seems worth while to develop the theory of the experi- 
ment in a general manner. We shall, however, retain the assumption 
that the conductivity of the material of the mass B is so good that the 
whole of the mass may be assumed to be of uniform temperature at 
each instant. 


DEVELOPMENT OF THE SOLUTION 


Convenient Relations 


It will be convenient to note the following relations between the 


constants of the system. 
We define \, as 


k 


(5) 


Xa 
Thus \, is the reciprocal of the “time constant’’ of the approximate 


solution, as exemplified by (4). 
We define R as 


— heat required to raise S one degree centigrade 
~ heat required to raise B one degree centigrade 


R= pol (area of contact) _ pol 
~ Ms (area of contact) Ms 
Again, since a? = k/po, we have 
_ kpo _: 
Msa? = Msk 


_ k Msl 
Msl pa 


Analytical Procedure 
Let us write 
= Ao + + Aox? + agx? + 


where the a’s are functions of the time. Then (2) demands that 


(do + Gyx + + + + ---) 
= (2a, + 3-2ayx + 4-3agx? + 5-4dagx? + ---), 


= 

3 

(7) 

Also 
k ral? 

a? = (9) 

po R 

10) 

(10) 

a 
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Since this must hold for all values of x, we require 


ay 
| 3a? 

a; 


5-402 


Qo = 


Thus, from (10) 


Gox? | doxt 
6 = + + + + 


a,x? a,x? 


+ ayx + | 3a? 


Now condition C requires that 


+5 +5 + =(. (13) 


We now try as a solution of this equation 


= Ae™*: a, = — Be-! (14) 


where A and B will become subsequently related to fit condition (3), 
or (3.1), whichever we choose to adopt. 


yielding 


(16) 


tio 
& 
d'*do 1 
aB 
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or, using (9) 


4| he (2 )") 0 
cos = A Rr sin = VU, 


= and A (>: ) =G (18) 


Writing 


ytany = 1/G. (19) 


B/A, and so G, become determined so as to fit condition (3), or (3.1), 
in accordance with our choice, and (19) then determines the possible 
values of ¥ and so, through (18), the possible values of A. G will, in 
general, involve \.. We shall return to these matters. For the moment 
we proceed to the general expression for @ as a function of x and 1. 


General Expression for 6 as a Function of x and t 
Writing \, for a typical value of \, and A, and B, and G, as the asso- 
ciated values of A and B and G, we have, from (14), 


G 


—hat 
(20) 


= Ae; @, = — Bert = — 


as the associated possible values of a) and a, satisfying (13), and so 
condition C. Moreover, any values of d») and a; given by (20) satisfy 
(13). We can, therefore, write (12), our solution for 6, in the form 


x 
0 a 


Using (9) and (18) 


l l 


where G, is to be determined by the condition (3) or (3.1), as the case 
may be. This gives@ = Oatt = ~ and sosatisfies (D). It remains to 


369 

2 

(17) 

6= = 

0 

a 

2 

6 
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determine the A’s so that they shall satisfy (B). Thus, we must have 


A, cos — Gy, sin = (22) 


We shall presently return to the problem of the determination of the 
A’s, but shall first discuss the forms of the G’s and the y,’s. 


FORMS OF THE G’s FOLLOWING FROM EQS. 3 AND 3.1, RESPECTIVELY 


We can proceed either by use of (21) or, more simply, by use of (12). 
Equation 3 requires that 
= ka, 


or, using (5) 


ao = 
Using (20) 
i, 
0 0 a 


Since this relation must hold for all values of ¢t, it must hold indepen- 
dently for each s. Hence 

G, = =. 23 

(23) 

Thus the boundary condition embodied in (3) leads to a value of G, 


which is independent of s. 
Considering now the requirement embodied in (3.1), we have, 


using (12) 


— — ka, = 0. 
Using (5) 
do — — d,/a, = 0. (24) 


Using (20), and observing that (24) must hold for all values of s, 
individually, we have 


RM 
kRu\ 
1 
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Observing from (18) that ¥, = (RA,/Aa)!?, Eq. 25 becomes 


1 
~ + ky2/Rhl) 


(26) 


G, 


so that now G, does depend upon s. 


THE NATURE OF THE j,’s AND THE },’s 


The y,’s are determined by (19) which yields, through (23), for the 
condition of Eq. 3, 


y,tany, = R (27) 


and through (26) for the condition of Eq. 3.1 


Rul) 


y,tany, = R (1 + (28) 


The X,’s are determined immediately in terms of the y,’s from (18). 
Considering both (27) and (28), it is obvious that the smallest y, 
that is, Yo, which we shall call 8», lies in the first quadrant, and succes- 
sive B's in the 3rd, 5th, 7th, etc., quadrants. The y,’s will therefore be 
of the form 


Yo= Wi 
with all the 6’s less than 7/2. 


¥2=2r+6.; sr +8, 
Thus (28) assumes the form 


k(sm + 8,)? 
Rhl 


(sr + B,) tanB, = R [ 4 (29) 


Case Corresponding to Equation of Condition (3)—-Absence of Thermal 
Resistance at the Junction of S and B 


For this case h is infinite, and (29) assumes the form 
(sm + B,) tan B, = R 


so that the 8, diminish with increasing s and, in fact, for large s, assume 
the form 


(30) 


R 
B, = tan B, = —. 
ST 


Thus, for large values of s, ¥, assumes the form 


= + R/s?x?*) 


or, to high approximation 


for large s. 
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= 

(31) 

= ST 
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From (18) and (31) 


+ R/s*n?)?, 


Rd,/a 


Thus, for large s, 


Ne. (33) 


Even for R as large as 0.5, and s = 1, 4, > 20\,. The result is that in 
our expression (21) for 6, only the term involving the smallest \,, which 
we shall call Ao, is significant in the sense that if, for example, e~t is 0.5, 
e~* will be less than (0.5)*°, and soon. Thus, unless the A,’s for large 
s grow to enormous sizes, which they do not, but, in fact, decrease with 
increasing s, we need concern ourselves only with the Ao, except very 
near tot = 0. 

It is now of interest to ascertain how near \o is to Aa. The value of 
Bo (corresponding to Wo) is given by 


Bo tan Bo = R. (34) 


For the case where R is small, it is necessary to recognize, using (18) that 

Bo( = Yo) = (35) 
where \» here refers to the smallest of the X’s. 
Ryo 1/2 
) 
For very small values of R, we can thus write 

yielding \o/A. = 1 for the limiting case where R = 0, as might be ex- 


pected. For the case where R is small but not zero, we may expand the 
tangent in (34) and write 


R = Bo (a + + ae (14 4 ). (36) 


The first approximation gives Bo? = R. Substituting this value for 
Bo? in the second parenthesis of (36) 


= 
2 
R. 


(RXo/da)!? 


Thus 


= 
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Hence, using (35) 


Bro) = R 


Xo 

1+ R/3 + R*/45 
If R ~ 0.2, Ao will be less than \, by about 7 per cent. The approxi- 
mate thermal conductivity is determined by writing 


do = k/Msl, 


where Xo is the experimentally determined value. The more correct 
expression to yield k& in terms of the experimentally measured Xp is 


k 

he = 38 

1+ R/3 + Msl(1 + R/3 + R2/45) 

Case Corresponding to Equation of Condition (3.1)—-Thermal Resistance 

at the Junction of S and B Recognized 

We must now use Eq. 28 which leads to (29). It no longer results 

that the 6’s necessarily decrease with increase of s. In fact, if # is very 

small, so that the thermal resistance at the junction is large, (29) 
assumes the form 


ti (sm + 
an B, = hl S17 B,) 


so that, with s infinite, tan 8, is infinite and 8B, = 1/2, which is the largest 
value of 8, attainable, since it must be in the first quadrant by definition. 

If, however, we consider the case where the magnitude of h is such 
as to permit Bo to be significantly less than unity, so that we may use 
the first two terms of the expansion of tan Bo, we shall have, from (28), 


(a+) = 6 (1 +4 = R + Bo. 


Our first approximation to 8, is given by 


1 — k/hl 


Substituting this for 8? in the second parenthesis, and neglecting terms 


: 
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involving products of R and k/hl, we readily find 


so that, from (35) 


1 

We see, as might be expected, that the effect of the thermal re- 
sistance at the junction is to diminish the difference between \ and Xo. 

Again, if \ is the experimentally measured value, the expansion for 
determining the conductivity k should be 


k k R 


In principle it is possible to determine 4 by performing two experi- 
ments with different values of /. In this way, two values of Ao are ex- 
perimentally measured, and provide us with two equations from which 
k and h may be obtained. 

While expansion of tan By has been invoked in the foregoing dis- 
cussions of both cases, that where thermal resistance is ignored and that 
where it is included, the use of such expansions is unnecessary. We can 
always handle the equations concerned graphically. We have used the 
expansions mainly to provide ease of description. 


THE DETERMINATION OF THE 4A,’s IN EQ. 21 


While the complete solution for @ as a function of x and ¢ is not 
necessary for the practical purpose of determining the relation which 
governs the measurement of k, it is of general interest to complete the 
solution expressed in (21) by outlining the processes which may be used 
for the determination of the A,’s. The fundamental equation for this 
purpose is (22). 

Of course, we can always adopt the procedure of starting by limiting 
the infinite sum on the left hand side of (22) to a finite sum of N terms, 
writing down N equations to provide a fit of the series to N points from 
x =0 tox =/. We then have N linear equations to determine the 
corresponding N coefficients A,, the results, of course, being only ap- 
proximations. By repeating the process with a larger number of terms, 
we can trace the progress of the approximation to any degree we desire ; 
but this process is laborious, and unimaginative for the purpose of 
bringing out the controlling influences of the various magnitudes in- 
volved. However, other methods of determining the A,’s present 
peculiar difficulties in the present instance. 


R 
= 
ki 63 
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If the expression on the left hand side of (22) were in the form of a 
Fourier series, with independent constants for the sine and cosine 
terms, the constants to be determined so as to produce the right hand 
side, all would be straightforward and simple as regards the determina- 
tion of the constants. However, in the first place the constants asso- 
ciated with the sine and cosine terms are bound by a common A,. 
Moreover, on account of the rather complicated nature of the ~’s 
occurring in the y,’s, the coefficients of the A,’s are not even ‘“‘normal 
functions.”’ For this reason, our determination of the A,’s must follow 
a line which is dictated by these circumstances. 

Recalling our discussion on the form of the G,’s, we infer that, for 
any boundary condition at x = 0, there is an infinite set of \’s, denoted 
by the \,’s, and an infinite set of G,’s related to the \,’s by the conditions 
which, in line with (19), are 

y.tany, = 1/G, (39) 
where, in line with (18) 
= (40) 


With this understanding, we can write (22) as 


A, cos 3 — cot y, sin y, | = Oy 


that is, 


us multiply by sin Y» (3 


sin 1) (7 1) dx = — sin vm (F 1) dx 


4 
i cos(¥m +W »s) 1) — cos(Ym — We) i) = [cos im (7 


. Ad sin(Ym + Ws) sin(¥m — ¥s) _ 
2siny, [ Vm + Vs Ym 7 Vs | Ym 
A,l [ Sin Ym COS Ws + COSYmSINY, SiN Ym COS Ps — COS Ym Sin 185 
25 = — (1 — cos Ym) 
2 sin Ym 


- 1) and integrate from x = 0 tox = 


Ym Vs Ym ~ Vs 


Using (39) 


> 
sin 
a: 
jin 
l 
Let 

0 
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For an assigned m, this expression holds with the sum extending to 
all A, from s = 0 to #, except that the contribution from A,, to the left 
hand side should be specially evaluated in view of a certain indeter- 
minancy in the term involving ¥, — y¥. in the denominator. To this 
end we omit the term involving A,, on the left hand side of (42) and 
replace it by J», where 


which is the appropriate contribution, as may be seen by (41). We 
thus have 


Ja = sin [cos (7 1) 1|ax 
~ 2sinynl Wea NJ “Jo” 2sinynl 2%n 


Equation 42 thus yields 


+ Gis — Gs 


sin Wn | 


Ym + Vs Ym 


Mn inva + = 
Gn SIN Wm (1 COS Wm). (43) 


2 


lA, + GwWs — 
* sin Wn 


Case Where Thermal Resistance at x = 0 is Neglected 
In this case, we observe that G, is independent of s and is equal to 
1/R. ‘Thus for this case, (43) becomes 


sin + OR (sin Vn + (1 — cos Ym) 


_ (1 — cos Hm Am 
A,= Vm ( sin Ym ) 2 + 


Again, using (39), this becomes 


_ 4% R ) 


eoMs 


R sin ¥ 
(44) 
R 
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Now referring to (41), let us apply the result at x = 0, in which 
case it yields 


(46) 


Thus, (45) yields 
266 


COS + R/sin? Pn)’ 


At this stage we may change the subscript m back to s, and write 


266 
~ cos ¥.(1 + R/sin? 


Recalling that y, = st + B,, we readily see that cos ¥, = (— 1)*8,, and 
sin? y, = sin? B,. Thus 


2(— 1)*6o 
A, = cos B,(1 + R/sin? B,)’ 


The substitution of (47) in (21), with G, replaced by 1/R for the case 
under consideration, completes the solution of our problem for the 
determination of @ as a function of x and ¢ in a form satisfying all the 
conditions A — EF. Referring to the development between Eqs. 39 
and 41, we see that, for the case under consideration, 


0 


Thus, since sin y, = sin (st + 8,) = (— 1)* sin 8B, 


sin B, cos B,(1 + R/sin? B,)’ (48) 


Comment Concerning the A,'s 

While the derivation of the A,’s as given in (47) seems perfectly 
logical, the procedure adopted is, as far as the writer is aware, rather 
unusual, and there would seem to be some mental satisfaction in apply- 
ing some independent test. Such a test would be to show that, with 


the A,’s as given by (47), ©&© A, = 00 as demanded by (46). The appli- 
0 


cation of this test seems to involve rather cumbersome procedures for 
the general case. However, it is possible to appy the test without much 
difficulty to the case where R is small. This we shall do. 


oA, = 6. = 
0 
m 

A, 
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Our basic Eq. 27 yields (30), which is 


(sr + B,) tan B, = R. 
For s = 0, this yields 
Botan Bo = R. 


For R small, both 8) and tan By are small, so that we may write 
3 2 
+ % = R. 


Substituting the first approximation 6)? = R in the second parenthesis, 
we have 


Bo? 
sin By = 


sin? Bo = 


R 


. = to the first power of R 
sin? Bo P 


R 
cos Bo = to the first power of R. 


Hence, from (47) 
265 


Ao = —R/2)(2 + 2R/3)" 


Retaining only first powers of R, we have 


R 
Ay (1 + 6 ). 


For s > 0, Eq. 49 yields as a first approximation 


R 


tan 8, = B, = 
Writing (49) as 


(sm + B,)B. (1 ) 


and substituting the first approximation for 8, in the two parentheses, 


| 
R 
1+ R73 
pe ( ) 
Bo | 3 = Bo 1 | 2 
R ( 
1 
(1+ R/3) is 
| 
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ST (1 + 1 + 3524? B, R, 


so that, retaining only the first power of R, we have 


B. = 


we have 


and to the same order sin? 8, = 8,2 = R?*/s?x?, and 
Ro _ 
sin?B, 
Again, 


cos 8, = 1 — 


Thus, from (47) we have, for s > 1, 
2(— 1)*6o 


| 
(1 — R2/2s?x?) ) 


A, = 


2( 1)*RO, 


R? ry 


Retaining only the first power of R 
_ 2(= 


A, 


So that, using (50) and (51) 
; 6 6 


It is not intended to imply that (46) holds for only small values of R. 
Our development requires, as a logical necessity, that it holds for all 
values of R. However, since it is not easy to verify this fact directly 
and independently, starting with (47), we have carried out this test 
for the algebraically simple cases where R is small. 


379 : 

1 - 

Ho. 

Thus 

A, = > 6, =_ lc IR ( A, =_— G5. (51) 

1 12 6 

ow 
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NOTE CONCERNING THE EFFECT OF THERMAL JUNCTION RESISTANCE 


The conductivity is calculated from \» which, to the approximation 
involved in retaining only the \» term in the solution for 6, at the junc- 
tion, gives 6 = @e~**'. The measured boundary temperature 6, differs 
from 6 by A@ cited in footnote (2). However, the relations in footnote 
(2) readily yield 


= /(1 = 6,/(1 hoMs/h). 


Thus the slopes of log 6% and log ®@ as functions of ¢ are the same, so that 
use of 6, instead of 6, does not affect the determination of \o. 


THE HUMAN OPERATOR AS A SERVO SYSTEM ELEMENT* 


BY 
DUANE T. McRUER! AND EZRA S. KRENDEL? 


PART I** 


SUMMARY 


This paper considers the role of human elements in certain closed loop control 
systems. A quantitative description of human dynamics useful to control system 
designers is essential to understanding and analysis of such systems. Accordingly, 
the human behavior description must be expressed in terms which are compatible with 
conventional descriptions of control system components. This compatibility is 
achieved by the use of a quasi-linear mathematical model for the human operator. 
The model is composed of two components—a describing function and remnant. The 
describing function, which for a linear system is identical with the conventional trans- 
fer function, is established to characterize that portion of the operator's output which 
is linearly correlated with his input. The input, upon which the describing function is 
based, is selected on the basis of a priori estimates of the nature of certain human non- 
linear behavior. Human output power which cannot be characterized by the opera- 
tion of the describing function on the input is designated as the remnant. 

After presenting the analytic basis for measurements of human dynamics, steady- 
state describing functions measured by the various experimenters in the field are dis- 
cussed and the adaptive, optimalizing behavior of the human operator is demonstrated. 
The remnants are also discussed and plausible sources for their origin are postulated. 

Knowledge of the range of parameter adjustment of which the human operator 
is capable in his adaptation as well as knowledge of his criteria for adjusting these 
parameters enables the designer to specify input functions and operator controlled 
dynamics compatible with both desirable human operator behavior and good system 
performance. By judiciously trading off system complexity against operator prefer- 
ences, while still making proper engineering use of the human operator's adaptability, 
a control system may be optimized for both performance and reliability. 


INTRODUCTION 


The control engineer achieves a logical connection between a system 
and its components by expressing component characteristics in terms of 
abstract models, such as transfer functions or computer setups, and then 
combining these models to obtain a description of system behavior. 
The use of such techniques in the design of automatic control systems 
has met with outstanding success. In this process, designers have 
tended to overlook one of the most versatile and, in many respects, most 


* This research was supported in whole or in part by the United States Air Force under 
Contracts AF 33(616)-3080, AF 33(038)-10420, AF 33(616)—2804, AF 33(616)-3610, and 
AF 33(616)-5822, monitored jointly by the Flight Control Laboratory and the Aero Medical 
Laboratory, Wright Air Development Center, Wright-Patterson Air Force Base, Ohio. 

1 Systems Technology, Inc., Inglewood, Calif. 

? The Franklin Institute Laboratories for Research and Development, Philadelphia, Pa. 

** Part IT will appear in this JouRNAL for June, 1959. 
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readily available and reliable components of all—the human controller. 
There are fairly broad areas where a suitable manual control system 
represents either the optimum solution or merits detailed consideration 
as a good possible choice. Before control engineering concepts can be 
applied to such cases, however, one must have a catalog of human static 
and dynamic characteristics, preferably expressed in terms readily 
understood by the control engineer. One of the principal purposes of 
this paper is to provide the reader with a synthesis, in control engineer- 
ing terminology, of the available data on the dynamic characteristics 
of human operators in certain continuous control tasks. This summary 
is to a large extent based upon McRuer and Krendel (1). Much of the 
available static characteristic data can be found in (2 and 3). 

Many agencies and individuals, both in the U.S.A. and abroad, have 
made substantial contributions to the state of knowledge of human re- 
sponse characteristics. In the past these data have been compared 
rarely if ever, since they were measured in different ways, different tasks 
were given the operator, diverse system inputs were utilized, and the 
data were presented in different forms. A review of the situation in 
the light of servo analysis indicates that consistencies exist for these 
data, from which useful generalizations can be hypothesized about the 
human operator dynamics. 

The human response measurements of concern here are those in- 
volved in closed loop control tasks which yield operator transfer charac- 
teristics for routine, essentially continuous functions such as stabiliza- 
tion and tracking. Even for these simple tasks we cannot expect to 
obtain a unique transfer characteristic since the operator’s responses 
will depend upon at least the following factors: 


1. The dynamic characteristics of the controlled elements, that is, 
the dynamics between the manipulated variable and the display. (The 
operator, to be successful, must adapt his transfer characteristics to a 
form required for adequate stability and performance of the entire 
closed loop system.) 

2. The type of input or forcing function driving the system. 

3. The individual reaction delays, thresholds, etc., of the human 
during the particular operation. (There are variations between indi- 
viduals and within individuals which are functions of both physiological 
differences and task experience. Any set of transfer characteristics ulti- 
mately proposed would have to allow for a variation in parameters of 
this type.) 

4. The motivation, attention, previous training, and general psycho- 
logical and physiological condition of the human at the time of the 
operation. 

These considerations specify a human mechanism which adapts itself 


’ The boldface numbers in parentheses refer to the references appended to Part I of this 
paper. 
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in some way as a function of the system inputs and controlled elements, 
and which may have wide individual variations in the values of the pa- 
rameters adopted. Individual parameter variations cannot be specified 
at present since there are too few data dealing with adequately large 
populations. There do exist adequate data to demonstrate the adapt- 
ability of human transfer characteristics to both forcing functions and 
controlled elements and such data provide a rationale for classifying 
past experiments. 

The selection of the form of the mathematical model to describe the 
operator’s transfer characteristic is based upon several practical con- 
siderations: the ease with which it lends itself to conventional servo 
analysis methods; the convenience with which experimental data may 
be obtained to define such a model; and the extent to which the measur- 
ing experiments are capable of generalization. These factors limit 
human operator models to the following : 


1. A linearized model consisting of an equivalent quasi-linear transfer 
function (which is correlated linearly with the system input or forcing 
function) and a remnant or residual output (which represents all output 
content which cannot be ascribed to a linear operation on the input). 

2. Simulation, using analog computer elements. 

3. A series of nonlinear cr piece-wise linear equations which can be 
solved by either analog or digital means. 


To date only the first two have been used successfully as the basis of 
experimental techniques. The observation that human operator con- 
trol responses are discontinuous has led to the conjecture that a non- 
linear relay servo model would provide a more natural characterization 
for human control behavior. Although certain preliminary experi- 
mental results indicated the existence of possible dual or multi-mode 
relay type servo behavior on the part of the human, other data shed 
doubt on this behavior (4, 5). The literature demonstrates that the 
model lacks convenient applicability to the analysis of continuous closed 
loop systems with continuous inputs (6, 7). 

To provide an adequate representation of the operator in a closed 
loop controller application, the operator mathematical models should 
provide answers to two questions: 

1. How does (will) the operator perform in a given situation, that is, 
what specific quasi-linear transfer function and remnant, or analog 
computer set-up, are required to describe his conduct ? 

2. How well can the operator perform and how do we make him 
behave that way, that is, what is the ‘‘desirable”’ operator characteristic 
and under what circumstances is it achieved ? 
The answer to the first of these questions is essential for manual control 
system analysis, while that of the second provides the starting point for 
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synthesis. The first question is answered by the formulation of an 
adaptive-optimalizing model developed to be consistent with all of the 
available data. This model can be used, with caution, for stability 
and limited performance predictions of manual control systems. The 
second question is answered to the extent of defining a desirable system 
form and the human operator’s response when operating with this form. 
This ‘desirable’ operator characteristic refers to a “best” characteristic 
within the adaptive properties of a so-called ‘‘average’’ human operator, 
and not to an ideal super-human-machine combination. What is de- 
sirable from the system designer’s lights may not be preferred from the 
controller’s viewpoint. 

In developing the answers to the above questions this paper is 
organized in four main parts. The first is primarily background infor- 
mation on the types of manual control systems of interest to this paper, 
and the general methods for the description of quasi-linear operation. 
The second presents a summary of presently available human dynamic 
response data, gathered from experiments directly involving the meas- 
urement of human describing functions from tracking systems which 
display the error alone. Past experimental efforts are reviewed, and 
the approximate mathematical models evolved to “‘fit’”” the data are 
presented. The third and fourth parts generalize from the experi- 
mental data. 

The hypothetical models discussed are in two categories: a general 
hypothetical model attempting a mathematical description of the human 
operator for various system considerations (that is, various controlled 
element configurations and types of system inputs) ; and a ‘‘desirable’’ 
human operator form determined by a minimum rms error criterion of 
system performance. Since the form of the highly adaptive human 
operator is determined by the control problem, the controlled elements 
of the system can be subjected to design change to elicit desirable human 
operator characteristics. By having both a general hypothetical oper- 
ator model and a preferred operator model, the control system analyst 
can evaluate the performances of particular system configurations with 
that of a more or less theoretical optimum, and can gain an idea of the 
type of modification required to the rest of the system if system error 
performance alone were the sole design objective. Since other factors 
such as reliability, maintainability and cost are often as important as 
dynamic performance, the designer must effect a compromise in order 
to achieve an over-all optimum. 


TYPES OF MANUAL CONTROL SYSTEMS AND METHODS OF OPERATOR DESCRIPTION 


Operator Environmental Factors Due to the Control System 


The fact of human adaptation makes description of the operator 
enormously complex when viewed in the large, and necessitates the 
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setup of simplified, constant situations within which one may have some 
hope of obtaining a simple behavioral model of engineering value. 

For the control situation to be analyzed with confidence the assigned 
task or goals of the human controller, and the constraints or rules govern- 
ing his interactions with the system, must remain reasonably constant 
over the duration of interest. In the present instance we are concerned 
with manual processes for minimizing visually perceived errors by exer- 
cising essentially continuous control so as to match visually presented 
input and output signals. This type of process is often called tracking. 

One can classify tracking problems for conceptual convenience. Our 
criteria for classification will be based on the type of information which 
the stimulus to be tracked and the controller's response present to the 
controller as a basis for future tracking decisions. In application, these 
criteria result in three limiting types of classifications: precognitive, 
pursuit and compensatory. These classifications arise as interactions 
among training, the geometry of the display, and the nature of the input. 
As such, they can readily exist in mixed forms. 


Precognitive. ‘This condition exists when the operator has complete 
information about the input’s future and a stimulus can trigger off a 
repertory of practiced, properly sequenced responses. In a sense, the 
situation itself may be the stimulus. Thus throwing a baseball at a 
target is precognitive behavior, as is steering a car out of a skid, or 
navigational flying under VFR conditions. The operator doesn’t need 
to maintain a frequent check on the individual responses in a sequence. 
Instead, end product responses are monitored; such as, “Did I hit the 
target ?’’, ‘‘Am I still skidding?” or ‘Have I passed over a given fix?”’ 
In that continuous close control zs not maintained on the perceived error, 
precognitive behavior is not tracking. However, we have included this 
classification because occasionally tracking approaches these conditions. 
One might characterize precognitive behavior as discrete control with 
long sampling intervals such that open loop control is in effect during 
these intervals. 

Pursuit. \n pursuit behavior past experience provides the tracker 
with information with which to predict the future input, but he must 
operate in a closed loop fashion with visual feedback about his responses. 
In this for.n of tracking the display effects of the operator's corrective 
responses can be distinguished from his input. These responses are 
displayed after they have been modified by the system dynamics. A 
gunner aiming for a moving target using an open sight is an example of 
pursuit tracking in which the system dynamics are those of the me- 
chanical system composed of gun, arm, and body. 

Compensatory. Compensatory tracking is the same as pursuit except 
that the visually displayed effects of the controller's responses are not 
distinguishable from the system’s input. Were our gunner aiming for 
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the target using a telescopic sight with a small field of view, we would 
have a compensatory situation. Since in compensatory tracking the 
visual display is the system forcing function minus the modified control 
response, the operator can determine the effects of control motion alone 
only under zero input conditions. 


The precognitive classification, representing a skilled operator re- 
sponse or calibrated, open loop system, is not readily amenable to servo 
analysis techniques. The pursuit system is a particular type of open 
cycle, closed cycle servo system, and requires a more complicated servo 
analysis (8). The compensatory classification, representing continuous 
close-control tasks, can be attacked by simple servo methods. The 
particular type of manual control system to be considered here falls into 
the compensatory category. 

The general nature of a compensatory display is represented in 
Fig. 1. The operator is presented with an input consisting only of an 


“TARGET” |, 
FOLLOWER’ 


Stationary Point) 


Error 


At) 


Fic. 1. Pure compensatory display. 


indicator showing the difference, or error, e(t), between a system forcing 
function, or command signal, 7(¢), and the system output r(t). The 
operator's task is to minimize the error signal presented by trying to 
keep the circle superimposed on the stationary dot. 

Besides the display, other fundamental factors in the control situa- 
tion are the type of element being controlled by the operator (such as 
an aircraft, automobile, etc.), and the actual means of exerting control 
(such as a control stick or wheel with their associated restraints, for 
example, springs, dampers, etc.). All of these characteristics will be 
lumped into the general classification of the “controlled element.” In 
many cases the controlled element can be described suitably by transfer 
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functions, though in others nonlinearities may be present and require 
description. 

The functional block diagrams of the simple control situation con- 
sidered, with the human operator as an element of a‘closed loop system, 
is shown in Fig. 2. 


SYSTEM 
FORCING OPERATOR OPERATOR SYSTEM 
FUNCTION ERROR VISUAL OUTPUT OUTPUT 

i(t) e(t) STIMULUS HUMAN c(t) CONTROLLED | 
OPERATOR ELEMENT 


Fic. 2. Functional block diagram of simple compensatory manual control system. 


Quasi-Linear Operator Description 

If the characteristics of the human operator for a given over-all task 
are assumed to be capable of quasi-linear description, the operator 
mathematical model will consist of a describing function plus an addi- 
tional quantity inserted as an input into the system by the operator. 
Then the functional block diagram of Fig. 2 can be made into an equiva- 
lent block diagram showing more detailed operations of the system. 
This is done in Fig. 3, where the symbols describe the various quantities 

Additional Quantity n.(t) | 


Inj b 
N. | 


| 
a | | 


OR WEIGHTING FUNCTION 
| c(t) oR C Ye(r) OR Y(jw) r(t) or R 


u(t) or! + E eV 
OR | 
FORCING “7” ERROR | Ww CONTROLLED 
FUNCTION Yp( jas) ELEMENT 


| TRANSFER FUNCTION | 


| HUMAN OPERATOR 


Fic. 3. Block diagram showing human operator in continuous control task. 


as time functions and as Fourier transforms. ‘To simplify the structure, 
the dynamic characteristics of the display are lumped with the controlled 
element, and the actual forcing function is modified, if necessary, to an 
equivalent one. The linear transfer characteristics of the human oper- 
ating on his presented input are described by the weighting function, 
y,(r), or the Fourier transform, Y,(jw). Since all of the operator's 


i 
J 
ay as 
357 
| 
| 
: 
: 


388 D. T. McRuer E. S. KrenpeEv (J. F. I. 


output is not described by the action of a linear transfer characteristic, 
an additional term, ,(t), is added at his linear output to complete the 
description of the total operator output c(t). This additional term is 
called the remnant. The location of 7,(t) at the operator’s output is 
arbitrary, and does not necessarily imply that such a quantity is 
physically inserted at that point. 

The block diagram of Fig. 3 illustrates clearly the servo system 
characteristics of manual control systems in general tracking tasks. 
This allows us to apply the whole body of servo theory in our attack 
upon human behavior in such manual control systems. As servo sys- 
tems with single inputs, the compensatory system is of extremely simple, 
single feedback form. As a closed loop system, the operator’s output 
and system error are given in terms of Fourier transforms and transfer 
functions by, 

jo) + Ne( je) 


Clie) = Ga) V.Ge) 


(1) 


SYSTEM 


INPUT ——»—J_ CHARACTERIZED BY: /—»—OUTPUT 


CHARACTERIZED BY: h(r) CHARACTERIZED BY: 
WEIGHTING FUNCTION 
| i(t) as a time function |. + (t) as a time function 
H (jw) 
2. I (jw) as the Fourier TRANSFER FUNCTION 2. R(jw) as the Fourier 
transform of i(t) transform of r(t) 


Fic. 4. Linear system representation. 


Types of Mathematical Models Useful in Describing Continuous Man- 
Machine System Behavior 


A review of some of the mathematical models which are useful to the 
controls engineer as means of describing physical devices leads directly 
to the concept of a quasi-linear describing function plus a remnant as 
basic tools for describing human dynamic characteristics. This concept 
allows the retention of the linear methods so successfully employed in 
servo analytical techniques. The remnant so developed is seen from 
Eq. 1 to be an element in the system error and output equations. In 
many manual control systems, it can be the dominant operator output 
quantity in the generation of tracking errors. 

With physical systems containing only elements which behave in a 
way describable by linear constant-coefficient differential equations, the 
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appropriate mathematical model is the element’s weighting or transfer 
function. If the system were represented by the block diagram of 
Fig. 4, where the weighting function, /(7r), is the time response of the 
system when an impulse function is applied at zero time, then the 
relationship between the response and the input would be given by the 
so-called superposition or convolution integral, 


ry = f = ver. (2) 
If the Fourier transform is used, the transformation of Eq. 2 becomes 


where R( jw), (jw), and H(jw) are the Fourier transforms of r(t), z(t), 
and h(r), respectively. Laplace transforms can, of course, be used in a 
similar way. 

The introduction of the weighting function, or its transform, the 
transfer function, allows us to describe the performance of a linear 
system element completely. If the convolution integral is not valid, 
and the behavior of the system is a function of the particular inputs and 
initial conditions the analyst must define system operation by input- 
response pairs. Fortunately many nonlinear systems of interest have 
specific input-specific response pairs which appear to be very similar to 
input-response pairs for linear systems. Thus the performance of some 
nonlinear elements, for particular inputs, can be approximated by a 
linear element called the describing function plus an additional quantity 
called the remnant.‘ 

A familiar example of this concept is the sinusoidal input describing 
function, which is derived from consideration of the harmonic response of 
the nonlinear element to a sinusoidal input at various frequencies and 
amplitudes. Consider a sine wave applied to the input of a nonlinear 
element having a single input and output. The output will very likely 
be a nonsinusoidal periodic wave with the same period as the input 
wave. If the output waveform is analyzed in terms of its Fourier com- 
ponents, the fundamental component will bear a relationship to the 
input sine wave which can be described in amplitude ratio and phase 
angle terms. The describing function will be the ratio of the funda- 
mental to the input in the same way as used in a linear system. The 
remnant will be composed of all the higher harmonics. 

Since the describing function concept is based upon a particular 
input type, describing functions can be defined for a particular non- 
linear element simply by considering different types of inputs. The 


* While the notion of molding a nonlinearity into an equivalent linear element is quite old 
and has been used by many writers, the first instance of major exploitation of the technique 
was probably by Kryloff and Bogoliuboff (9). 


/ 
; 
3 
: 
q 
: 
5 
if 
5 


390 D. T. McRuer anv E. S. KRENDEL (J. F. I. 


sinusoidal input describing function already mentioned is the most 
widely known example. Where the type of transient input is known 
and simply described in analytic form, a transient describing function 
may be used; for example, the step input describing function. In 
cases where system input forms may not be simple steps or sine waves, 
a more useful describing function may be one based upon statistical 
inputs (10,11). Such a describing function, particularly for inputs 
characterized by Gaussian amplitude distributions, is useful in nonlinear 
control problems where the inputs expected in the application of the 
control system are statistical, and can be considered as a basis for a 
fairly general type of equivalent linearization. The Gaussian input de- 
scribing function will be shown here to be of particular value in describ- 
ing human response. 

As a general comment on describing functions, it should be noted 
that, since the inputs are different, the only thing the various describing 
functions have in common is the nonlinear element they intend to ap- 
proximate. As the elements become less nonlinear, the various describ- 
ing functions for a given element tend to approach one another. When 
the nonlinearities are entirely removed, all of the describing functions 
become identical to the element’s transfer function. 

Alternate approaches to the nonlinear problem exist. The most 
practical alternate for human response measurements is time domain 
analysis by the use of an operator analog comprised of analog computer 
elements. Using the very wide range of nonlinear and linear analog 
elements available, a computer setup can be made by cut-and-try pro- 
cedures and adjusted until the ‘‘analog operator’ responses to particular 
inputs are similar to those of the actual operator. If an all-compre- 
hending analog model were ever achieved, the input type would not be 
important. To achieve a practical and reasonably simple computer 
setup in the case of the human operator, however, a separate analog 
is required for different inputs. The analog model technique is of great 
value in instances where a point by point prediction of operator response 
is desired, for studies including nonlinear control effects, and as a means 
of providing insight into some of the types of nonlinear behavior which 
might occur in the operator. When nonlinear or time-varying transfer 
characteristics become exceedingly important it is probably the only 
practical approach. The analog technique is fairly straightforward, so 
no further discussion of the general method need be included here. 


Quantities Involved in Human Response Measurements 


Human operator characteristics are strong functions of the type of 
system input, so the describing function concept is necessary if we wish 
toretain and apply linear methods of analysis to manual control systems. 
Almost all experimental research to specify operator dynamics in manual 
control tasks has been based on this approach. A key factor in the 
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determination of the operator’s response characteristics is the general 
predictability of the system input. Specifically, the operator’s mode of 
response with either a unity controlled element or a thoroughly learned 
control element is different for single sine waves (which are essentially 
completely predictable once recognized), step functions (which have 
known final values immediately after initiation of the stimulus), and 
“random appearing functions.”” The random appearing inputs are not 
restricted to functions which are known only by their statistics, but 
include those which can be made up by as few as three nonharmonically 
related sine waves. In most of the experimental work, the random 
appearing inputs have been either such a sum of sine waves, or time 
functions which had essentially Gaussian amplitude distributions. 
Therefore the types of describing functions of general interest in human 
response measurement are those for step, sinusoidal, and Gaussian in- 
puts. All of these describing functions fall into a category which can 
be called quasi-linear because they tend to be linear under a fixed set 
of conditions, such as given inputs, yet nonlinear when changes in these 
conditions are considered. Since the statistical forms are usually the 
closest to reality in practical control systems, we shall confine our 
attention to transfer characteristics derived for random appearing 
inputs. When random inputs are used, the most common techniques 
of describing the characteristics of the time stationary inputs and out- 
puts are by means of correlation functions or power spectral densities. 
The describing functions are obtained from cross-correlation or cross- 
spectral density functions, and the remnant and system forcing functions 
are normally expressed as spectral densities. While most of the under- 
lying theory on the measurement and interpretation of these quan- 
tities is beyond our present scope (see, for example, 12 for details), 
it is necessary to develop some relationship among the fundamental 
quantities which are used to measure and define the operator's response 
characteristics. 

The various quantities in the compensatory system of Fig. 3 which 
are important in the determination of operator characteristics are given 
as a nomenclature summary, in the accompanying tabulation. 

The operator’s output and the system error for the compensatory 
system were previously given in Eq. 1. The cross-spectral density 
functions ®,, and ®,, are defined as follows (where the asterisk indicates 
a conjugate), 


= lim = lim jw)C(jw) (4) 
1 Tox 
= lim —[I*E] = lim a. L1( jw) E (jw) |. (9) 
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Expressed as Expressed in 
a Time Expressed as a Spectral 
Quantity Function Fourier Transform* orm** 


System forcing function i(t) I(jw) or I 


Error signal of the system ; e(t) E(jw) or E 
input to the operator 


Operator output c(t) C(jw) or C 
System output response r(t) R(jw) or R 


Remnant at the operator’s n.(t) N.(jw) or 
output 


Operator’s weighting and Y p(t) Y,(jw) or Y, 
describing functions 


Controlled element weighting ye(r) Y.(jw) or Y, 
and transfer functions 


Closed loop weighting and de- h(r) IT(jw) or H 
scribing function between 
forcing function and opera- 
tor’s output 


Cross spectral density of fore- 
ing function and operator 
output 


Cross spectral density of sys- 
tem error and operator 
output 


Remnant expressed as a n,'(t) N.’(jw) or Pan 
closed loop quantity 


* So that the various time functions, which may be stochastic, may be considered to have 
Fourier transforms, it should be understood that the time function used in the transform is 
identical to the actual time function in the intervals — T <t < T, and is zero elsewhere. 
The Fourier transform of a typical signal, f,(t), will then be defined as 


F(jw) = where f(t) =fil); —T<t<7 
f(t) = 0 elsewhere. 
** The power spectral density of a typical signal f,(t), is given by 
dy, = lim = F(jw) F( — jw). 


The cross spectrum of two signals, f;(¢) and g,(t) will be 


= lim = F( — jw)G (jw). 
1 


For more details on these quantities see, for example, 13. 
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Using Eqs. 1 and 4, the cross-spectral density ®;, is seen to be 


Ihe quantity, lim T [I*I ] is recognized as the power spectral density of 


the forcing function, The quantity lim [N.I*] is similarly 


T+ 0 


recognized as the cross spectral density between the forcing function and 
the remnant, ®;,,(jw). This cross spectrum is zero, since by definition 
the remnant has no linear coherence with the forcing function. Then 
Eq. 6 becomes 


(7) 


Similarly, it can be shown that the cross spectral density between the 
error and the input is 


1 
[1+ Y,¥.] 
The open loop describing function Y, can be obtained by dividing 
Eq. 7 by Eq. 8, and is given by 


(8) 


= —. 


(9) 


Equation 9 is fundamental in the measurement of the transfer charac- 
teristic of stable system elements. It is applicable to either random or 
periodic time functions, which must be stationary. If the system is 
linear and the forcing function a sine wave, the determination of Y, re- 
duces to a conventional frequency response measurement. 

The spectral density of NV, is often an important ‘“‘property”’ of an 
equivalent linear system. To find this quantity we can form the oper- 
ator output power spectral density, noting, as before, that the forcing 
function remnant cross spectra are zero, 


lim 3 [C*C] = lim = 
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In most cases it is easiest to consider the total output power in terms of 
that portion linearly coherent with the forcing function (the first term 
above) and a remainder or ‘“‘closed loop” remnant. Since Y,/(1+ Y,Y.) 
is recognized as a closed loop describing function, H7, and defining 
+ Y,Y.) |? ®an, as the closed loop remnant power spectrum ®,,, 


we can write Eq. 10 as 


The possible sources and points of injection of the actual quantities 
composing the closed loop remnant will be discussed later. For the 
present, the term ®,,, is selected to represent, without regard to origin, 
that portion of the operator's output power density which is not 
linearly coherent with the forcing function. 

Besides the describing function and the remnant, it is desirable to 
have some alternate means of relating the amount of operator output 
power “‘explained”’ by the action of the quasi-linear describing function 
to the total output power. We can define such a quantity by forming 
the ratio of the linearly correlated output power and total output power 
spectra. This fraction is the square of the linear correlation, p, and is 
given by 


LT \ 
(12) 


The linear correlation can also be expressed as 


& 


It is apparent from Eq. 12 that for a linear system (remnant, ®,,, equal 
to zero) the value of p would be unity. In general, the value of p isa 
good measure of the importance of the remnant. As will be seen 
later, the linear correlation is not necessarily a good measure of linearity, 
although a value near unity does imply that the system is ‘“‘almost’”’ 
linear. 

All of the quantities presented above are applicable to most general 
systems where servo analysis methods have a reasonable chance of 
success. They are essential, of course, in attempting to describe human 
operation because of the enormous adaptability of man as a control 
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SUMMARY OF HUMAN RESPONSE DATA FOR COMPENSATORY TASKS 


Introduction 


A. Tustin first noted that operators in manual control systems with 
random appearing forcing functions exhibit a type of behavior directly 
analogous to that of the equalizing elements inserted into servo systems 
to improve the over-all performance of the system. By its nature, this 
equalization appears to be a more or less continuous operation per- 
formed on data observed by the external senses. Since the operator 
must adapt his equalizing behavior to cope with controlled elements 
having widely diverse dynamic characteristics, it is reasonable to classify 
describing function experiments on the basis of the controlled elements 
which ‘‘mold” the operator’s characteristics. The discussion of oper- 
ator describing functions in compensatory tasks has therefore been 
divided into two major parts: the first will treat past experimental 
describing function data on compensatory systems with unity gain 
controlled element dynamics, that is, Y. = 1, and the second will con- 
sider describing functions for operators in systems with more complex 
controlled element dynamic characteristics. To complete the charac- 
terization of the operator, data on the remnant are also presented, 
though pertinent data are sparse. 


Operator Describing Functions in Compensatory Tasks with No Controlled 
Element Dynamics 

The experimental data considered for Y, = 1 were taken directly 
from the efforts of three investigators and their colleagues, specifically 
J. I. Elkind (14, 15), L. Russell (16) and E. S. Krendel (17). The pri- 
mary experimental variables in these researches were the type of con- 
troller and the forcing function characteristics. The details of the 
forcing functions used are summarized later in Tables I] and III. These 
experimental data have been replotted in Bode diagram form (when the 
original data were not in this form) and describing functions for the 
operator have been obtained from curve fits to these plots. The 
describing functions as derived from these curves are the basis for 
making analytical estimates from which one can determine the influence 
of the forcing function upon the human operator’s dynamic charac- 
teristics. To the extent that the curve fitted analytical model approxi- 
mates the operator, we have insight into the laws which actually 
govern the over-all behavior of the operator. 

The most desirable types of linear models are transfer functions 
made up of ratios of rational polynomials, with a pure time delay term 
also allowed. Such transfer function forms are favored because they 
are simple, well understood, and adequate for approximating any of a 
large number of frequency response characteristics. Unfortunately, 
since the experimental data covered a restricted frequency range, with 
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poor resolution, it is possible to fit functions to the limited number of 
data points to an arbitrarily defined degree of accuracy. A painstaking 
and accurate fitting procedure is not warranted because of both the 
measurement and data reduction errors and the variation between sub- 
jects whose data were averaged. The appropriate model to use in 
fitting a set of data is, therefore, the simplest transfer function form 
which is reasonably consistent with all of the data trends and is ap- 
propriate for both amplitude and phase characteristics. The need for 
consistency between amplitude and phase forces compromises in the 
estimates of parameters. 

As general criteria for the types of transfer function forms acceptable 
we list the following which are consistent with all known human operator 
data in compensatory tasks of various kinds. 


1. The operator’s describing function must approach zero at infinite 
frequency, that is, the transfer characteristics are fundamentally those 
of a low pass filter. 

2. The describing function must be finite at zero frequency because 
of physical limitations on the force-displacement capabilities of the 
operator. 

3. A pure time delay, e~**, must be included to account for a reaction 
time delay. 

4. When the experimental closed loop system is stable, the approxi- 
mate fitted describing function must define a stable closed loop system. 


As an example of this fitting procedure, we shall consider the case 
of Elkind’s data for forcing functions with rectangular spectra. (The 
final describing functions for this case are summarized in Table III.)® 
The simplest transfer function form that satisfied the foregoing require- 
ments for these data is 

Ke-nie 
Y,, = (14) 
 (Tojw + 1) 
where K is the d-c. gain, 7; a pure time delay, and 7) is a time constant 
representing a low frequency lag. Since the pure time delay e~"* does 
not affect the amplitude characteristics, the amplitude ratio curves are 


fitted by 


Approximate values for K and 7) can be determined from the raw data 


5 The original data upon which all of our discussion is based can be found in the original 
sources in various forms. They are also presented in (1) in Bode diagram form. Because of 
the enormous amount of data involved (Tables II and III summarize over forty separate ex- 
perimental averages derived from over 300 individual runs), it is patently impossible to show 


even a small variety here. 
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amplitude ratio plots with the chief criterion the visual appearance of a 
good fit. 

The phase associated with the term (7)jw + 1)~! can then be sub- 
tracted from the raw data phase to isolate the phase lag due to e~"*. 
On a linear plot of phase versus frequency, this residual phase should be 
approximately linear with a slope of — 7;.. When the residual phase is 
not approximately linear, it is necessary to select other values of 75 
(and hence K), and repeat the process until satisfactory results are 
attained. 

While the low frequency approximation of Eq. 14 fits much of 
Elkind’s data quite well, an analytical closed loop transfer function 
developed from this equation, with numbers derived from the experi- 
mental data points, would indicate an unstable condition at frequencies 
well beyond the measurement bandwidth. This was the case for all 
of Elkind’s lower cutoff forcing functions. Since there was no evidence 
of a closed loop instability in the experimental data, the operator pre- 
sumably introduces more attenuation at high frequencies than is indi- 
cated by Eq. 15. Although the experimental data do not define the 
high frequency portion of an approximate transfer function, it is ap- 
propriate to “‘stabilize’’ the closed loop system by adding a lag to Eq. 14 
and modifying the value of 7; to 7, that is, changing Eq. 14 to, 

(Tojw + 1)(Tijw + 1) 

If 1/7, is much greater than any of the frequencies measured, then the 
low frequency behavior of the transfer function of Eq. 16 can be made 
the same as that of Eq. 14 if 7, + 7 is approximately equal to 7;. This 
addition of a simple lag, 7,jw + 1, is not a unique way to place a 
stability bound upon the possible higher frequency behaviors; any 
combination of lags and leads leading to a stable system would be suit- 
able if it satisfied the relationships, 


Y,(jw) = (16) 


T+ ~7;- = Ti > (17) 
i k Ti 
where the 7; are lags, the 7), are leads and w,, is the forcing function 
cutoff frequency (or an approximation thereto) in rad/second. ‘The 
single lag represented by 7, is, however, the simplest form which can 
be added to Eq. 14 to assure stability and consistency with the measured 
low frequency response. Even with this procedure it is unfortunately 
not possible to find unique values of 7 and 7); from the low frequency ex- 
perimental data. 
The results of processing the various experimental data for systems 
with unity controlled element dynamics in the fashion described above, 
are shown in Tables II and III.* These tables delineate the experi- 


6 The symbol s has been used in all summary tables for jw. 


a 
| 
bis 
- 


398 D. T. McRuer anv E. S. KrenpDEL (J. F. 


mental conditions, the frequency ranges of measurement, the average 
linear correlations, and the approximate describing function fits to the 
data. Table I is also included to provide a summary of operator re- 


Caste 1.—Summary of Operator Describing Functions in Compensatory Tasks. 
Simple Tracker, Y. = 1, Simple Forcing Functions. 


“Best Fit’’ Human Operator Investigators 
General Control Task Transfer Function and Remarks 


Simple following with pencil, Closed loop transfer function Goodyear (35) 
stick, or wheel (forcing eT eT also: 
function a) 2¢7Ts? Tvs+1 Searle & Taylor 
r=0.5 K +6.8-7.8 Cheatham 
=0.25 sec. T =0.042 sec. Ty =1/K Mayne 


Ellson & Hill 


Simple following with pencil “‘ Nonsynchronous” response, Goodyear (19, 35) 
Mayne, Ellson & 
Gray, Cheatham 


or wheel (forcing func- closed loop 


tion b) 1 
T2538 2Ts 
K + K °K +1 


with constants as above 
“Synchronous” response, 


closed loop 


K (w) 


Simple following with pencil, Same as for Goodyear (19, 35) 
stick, or wheel (forcing single steps Mayne 
function c) Though 7 can approach 0.20 sec. as small range effect ; 

steps get closer together and Searle & Taylor 
appear more random Ellson & Wheeler 
Slack 


Types of forcing functions: 


_ 


a. Step functions b. Sine waves and square c. Sequence of steps 


waves 


sponse to simple forcing functions. These are not considered in detail 
here, but have some general interest and are also of some value in ob- 
taining our hypothetical human behavior models. 

Tables I] and III shows a wide variety of forcing function conditions 
represented by their power spectral densities, ®,;. The broadest 
coverage is due to Elkind (Table III), who had remarkable success 
establishing relationships between the fitted describing function param- 
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These are also given in Table IIT. 
ships. While these attempts were not ¢ 


justment to minimize the rms error (13) 


TABLE III.—Summary of Operator Describing Functions 
Forcing Functions. Simple Tracker, Y, 1, Random A 
Taken from Elkind, in a Control Task Consisting of 


“Best Fit'’ Human Operator 
Transfer Function 


Ke-™ 
(Ts +1)(Tws+1) 


Type of Forcing 
Function 


Function made up 
of 40-120 sinus- 
oids, giving any 
desirable rectan- 
gular spectra 1’ 
rms amplitude 
(see g below) 


where: 
K 
Actual fitted data: 


=2.2/ feo? ; 


Forc- 
ing 
Funce- 
tion 


| R.16 


K (db) 
34.5 


| R24 315 104 314 6.22 
R.40 22.5 .785 12.3 
R.64 50 150 1.73 30.3 
R.96 6.5 139 3.65 x 
| R16 -06 122 3.77 « 
| R24 —3.0 116 1,885 «x 
Approximate white ke 
noise through 1st, | T 
2nd, and 3rd order | (Tis+1) 
binomial filters as | where: 
shown below. 3 
Highest fre- K +=+— 
quency in forcing 
functions was lor 
2.88 cps. 1” rms 
amplitude (see h | K+9.42T, 
below) of 
4 


Actual fitted data: 


Forcing 
Function 


and Filter K(db) 1/7, T 
Fl; 1st order 10.5 1.13 .139 
F2; 2nd order 25.0 0.314 126 
F3; 3rd order 33.0  .1885 .178 
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eters and other parameters defining the forcing function characteristics. 
Some work has been done (1) in an 
en/leavor to find underlying criteria which may lead to these relation- 
completely successful, one can 
show that the conventional servo synthesis method of parameter ad- 
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leads to analogous results. 


in Compensatory Tasks Showing Effects of 
ppearing Forcing Function. These Data 
Simple Following with a Pip Trapper. 


| 
Frequency Range of Human | Average Linear 
Operator Measurements | Correlation, pa 


Highest frequency equal | —R.16-.995 
to forcing function | R.24-99 
bandwidth, lowestfre- | R.40—-.995 
quencies given below R.64-.98 
in rad/sec. R.96-.92 

R1.6-.75 
R2.4-.58 


R2.4—.5 
Range Filter Pa 
Filter (rad/sec) Ist order 0.9 
Ist order .25 to 18 2nd order 0.96 
2nd order .25 to 12 3rd order 0.97 
3rd order .25 to 6 


R.16-.15 
R.24—.15 
R.40-.25 
T 1/Ty R.64-—.25 
R.96-.25 
| 
| 
; 
| 
| 
| | 
— = | 


May, 1959.] HuMAN OPERATORS IN SERVO SySTEMS 


TABLE II1.—Continued. 


Type of Forcing | “Best Fit’’ Human Operator Frequency Range of Human Average Linear 
Function Transfer Function Operator Measurements Correlation, pa 
24 to 144 sinusoidal me | 
components giv- (T1s+1) (Tws+1) | 
ing forcing func- 
tion spectra Jwhere Y(s) was translated to 
shown (No. of | the origin for B3, B8, B9, B10 B3 yo — -9 to | B3 oe 
components B4 6 ‘on 
_ | where: .6 to 8.7 B4 90 
parentheses.) 1” | K29.427,; K2—S_ to8 | Br 
rms (see 7 below) ae! BS 3to6 | Bg 07 
Actual fitted data: | BO 6to9 | Bo ot 
| B10 9to12 | B10 70 


Fore- | 
ing | 

Func- 
tion AK(db) 1/7, 1/Ty 
Bi 9.0 4.78 «© 153 | | 
B2 1.5 5303 107 | 
B3 —1.0 12.6 L .278 | | 
B4 —0.3 12.6 | 
B5 11.1 1.88 128 | 
Bo 17.7 1.00 .149 | | 
B7 23.2 88 17.8 .100 
B8 3.14 6:28. | 
B9 0.4 12.6 390 
B10 —11.8 2.82 


~ ESL 


88 


feo = .16, .24, .40, .64, 
.96, 1.6, 2.4 cps 


or 


Wee = 1, 1.5, 2.5, 4, 6, 
10, 15 rad/sec. 


(To be continued) 
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Simple Electronic Circuit Simulates 
Living Nerve Cell. simple elec- 
tronic circuit that simulates some 
functions of the individual biological 
nerve cell, or neuron, has been devel- 
oped at Bell Telephone Laboratories. 
Numbers of these artificial cells are 
being combined into experimental net- 
works that are roughly analogous to 
the nerve systems of the eye and ear. 

Bell Labs scientists are especially in- 
terested in discovering how visual and 
auditory nerves function and how their 
signals are interpreted by the brain. 
IKXnowledge derived from such basic 
research may lead, as it has in the 
past, to better and more economical 
communication. 

The gross function of the nerve cell 

transmission of electrical pulses in 
response to stimuli, and only to those 
stimuli that meet certain conditions 
has been simulated by a simple circuit. 
L.. D. Harmon, who initiated the proj- 
ect, has combined the cells into groups 
that simulate simple functions of the 
eye, and similar experiments with ear 
models are being started by W. A. 
Van Bergeijk. 

The neuron circuit, designed by R. 
M. Wolfe, fires electrical pulses of 
standard amplitude and duration, just 
as a biological cell usually does. If 
the circuit is driven by a constant 
stimulus, simulating receptor cells as 
in the eye or ear, trains of pulses are 
emitted. A higher intensity of excita- 
tion increases the frequency of pulsing ; 
and when the neuron is excited con- 
tinuously, the frequency of pulses can 
be made to decrease with time, exhib- 
iting accommodation as a living nerve 
cell does. 

Input excitation must, as in a bio- 
logical cell, surpass a threshold value, 
and the cell will integrate two or more 
input pulses below threshold value to 
cause firing. A particular input con- 
nection can also, while energized, in- 
hibit firing of the neuron by other in- 
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puts. Similarly, immediately after 
firing, the electronic neuron’s thresh- 
old rises to infinity and for a few milli- 
seconds no input signal can fire the 
neuron again. 

The circuit includes four transis- 
tors, thirteen resistors, and two capac- 
itors, mounted on a_three-by-four- 
inch, printed-circuit card for ease of 
handling. The pulse length it de- 
livers, about six milliseconds, is con- 
siderably longer than that of a bio- 
logical nerve, but can be shortened if 
desired. The cell has an integrating 
time-constant of two milliseconds and 
a refractory time-constant of about 
ten milliseconds, approximating time- 
constants of the biological neuron. 
Because the electronic inputs and out- 
puts are compatible, the cells can be 
assembled into chains and networks. 

Electronic neurons can be com- 
bined with photo-resistive cells to sim- 
ulate simple functions of nerves in the 
retina. Some receptors, known as 
“on” receptors, fire only when the 
light intensity they receive is increas- 
ing; “off” receptors fire only when the 
light is decreasing; and “during”’ re- 
ceptors fire while they receive a steady 
light. Flicker-fusion phenomenahave 
also been produced. In the human 
eye, these can cause a sequence of 
flashes to be seen as continuous illumi- 
nation; this property of vision is ex- 
ploited in motion pictures and TV. 

The mutual inhibition of cells in an 
array has been demonstrated experi- 
mentally. Some animals have been 
observed to possess this arrangement, 
in which a cell receiving a greater light 
intensity inhibits the firing of nearby 
cells that receive less light. This re- 
sults in local sharpening of image 
boundary detail. 

It is hoped that further imitation 
of nerve functions, so far as they are 
known, and other experiments may 
lead to better understanding and pre- 
diction of neurological behavior. 
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A UNIFIED APPROACH TO THE THEORY OF 
SAMPLING SYSTEMS* 


BY 
R. E. KALMAN! AND J. E. BERTRAM? 


ABSTRACT 


This paper presents a new method for the analysis of sampling (or sampled-data) 
systems. Unlike methods currently in use, it is not necessary to assume that the 
sampling operations are synchronous, performed at constant rate, and representable 
by means of an impulse modulator. Several different types of sampling operations are 
considered in detail, with the analysis proceeding similarly in each case. The paper 
concludes with a brief study of the stability of sampling systems and a generalization 


of Floquet’s theorem. 


INTRODUCTION 


A dynamic system in which one or more variables are allowed to 
to change only at discrete instants of time (called sampling instants) is 
known as a sampled-data or more simply as a sampling system (1).° 
Sampling systems may be regarded as a subclass of nonstationary (time- 
varying) dynamic systems. Such systems arise naturally when a digital 
computer, a pulse or counting circuit, or any other discontinuously act- 
ing element is incorporated in a system. 


Conventional Sampling 

At the present time, most of the literature (1-5) dealing with the 
analysis and design of sampling systems is based on the following sim- 
plifying assumptions: 


(a) All sampling operations can be adequately described by intro- 
ducing a fictitious idealized component, called an zmpulse modulator (2) 
into the system. 

(b) All sampling operations are performed synchronously. 

(c) The time interval between successive samples is fixed. 


In addition, it is almost always assumed that the system is linear. 

A sampling system whose operation (within the accuracies demanded 
by engineering considerations) satisfies these simplifying assumptions 
will be referred to in this paper as a conventional sampling system. 

The analysis and synthesis of conventional sampling systems have 


* This work was done while both authors were in the Department of Electrical Engineering 
at Columbia University. The work was supported in part by the IBM Research Laboratories. 

1RIAS, Baltimore, Md. 

21BM Research Laboratory, Ossining, N. Y. 

3 The boldface numbers in parentheses refer to the references appended to this paper. 
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been well developed and standardized in recent years with the aid of the 
z-transform method (1, 3-5). 

Recently there has been increased interest in systems in which the 
sampling operation does not satisfy one or more of Assumptions (a), 
(b), (c). Examples of such systems are: 


Nonsynchronous Sampling 

All sampling operations are repeated at the same rate but occur at 
different instants of time. Such problems can be analyzed with the aid 
of the ‘‘modified”’ z-transform method (1, 4, 5) due to Barker. In all 
except the simplest cases, however, the algebraic details of the analysis 
become exceedingly laborious. 


Multiple-Order Sampling 

All sampling operations in the system are performed synchronously 
and a group of sampling operations is repeated every T seconds, the 
intervals between successive samples in the group being unequal. For 
example, m**-order sampling means that once very T seconds samples 
are taken at times 
t=kT, RUE +1m-1, (RH+1)7,--- =0,1,---) 
where the 7; are arbitrary constants subject to the condition 0 < 7, 
<T, < +++ < ty < T. Such sampling operations may arise when a 
digital computer is time-shared by several feedback control systems. 


Systems with multiple-order sampling have not been studied previously 
as far as the authors are aware. (Added in proof: See 20, 21.) 


Multi-Rate Sampling 

There may be several sampling operations in the system with fixed 
but unequal sampling intervals. Kranc (6) has made an extensive 
study of multi-rate systems in which the ratios of the various sampling 
periods are rational numbers. His analysis is based on a special trans- 
form method. 


Noninstantaneous Sampling 

Here the sampling interval is too long with respect to the time-scale 
of the dynamic components of the system to be adequately described by 
impulse modulation. An approximate solution to this problem has been 
obtained by Kranc (7), using transform methods. An exact solution 
has been published by Farmanfarma (8). See also the writers’ dis- 
cussion of Farmanfarma’s paper. 


Random Sampling 
In some cases, the intervals between successive samples may be 
thought of as being selected according to some random scheme. Ap- 
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proximately random sampling may be caused by the desire to reduce 
susceptibility of the system to jamming or similar interference; by 
random time delays in the sampling operation depending on a physical 
measurement ; by a random time delay which must elapse before a digital 
computer again becomes available to carry out control computations 
required by a particular feedback loop of a multi-loop system; and 
finally by inaccuracies in equipment which is designed to sample at 
constant intervals. 

A comprehensive study of systems with random sampling has been 
carried out by one of the writers using mainly the methods of this paper ; 
these results are discussed elsewhere (18). (Added in proof : See 22.) 

Evidently, all five types of sampling mentioned above may be en- 
countered in a single system. 


The New Method 


The objective of this paper is to present a uniform method of 
analysis which can be applied in a straightforward fashion to sampling 
systems of the types mentioned. The same method has been used 
recently by the writers in connection with: 


(i) design of conventional sampling systems for ‘deadbeat’ re- 
sponse (9) ; 
(ii) design of an optimal relay servo operating with conventional 
sampling (10) ; 
(iii) design of conventional sampling systems to be optimal in the 
Wiener sense (11); 
(iv) systems with random sampling (18) ; 
(v) analysis of effect of round-off errors on the dynamic perform- 
ance of sampling systems incorporating a digital element (12) ; 
(vi) analysis of conventional sampling systems; and 
(vii) design of Wiener filters operating on sampled data (19). 


Transform theory as used in the past to study conventional sampling 
systems has proved to be either unwieldy or has required elaborate 
modifications before it could be applied to systems of the type men- 
tioned. By using the concepts of state and transition matrix (see the 
next two sections) it is possible to handle systems of the general type in 
a clear and uniform way. The new method yields simplifications even 
in the analysis and synthesis of conventional sampling systems. The 
method automatically eliminates one of the chief difficulties of the trans- 
form method, namely that with the latter it is difficult or cumbersome 
to obtain information about the behavior of the system at any time 
other than the sampling instants. The method is in accordance with 
modern trends in engineering science ; it devotes paper-and-pencil stud- 
ies only to the analytic aspects of system problems, leaving the drudgery 
of numerical computation to be performed by a digital computer. 


4 
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Finally, the formulation of system equations according to the point of 
view of this paper is a natural and necessary step preliminary to the 
study of nonlinear sampling systems. 


ELEMENTS OF LINEAR SAMPLING SYSTEMS 


Because of the generality of the method of analysis to be presented, 
careful attention must be paid to the definition of the individual ele- 
ments of the system. The elements of a sampling system to which the 
method applies are of the following type: 


Dynamic Elements 


These are elements whose present output depends not only on the 
present value of the input but also on past values of the input. We 
will be concerned with two particular classes of dynamic elements: 

Linear Continuous Dynamic Elements. ‘These are assumed to be de- 
scribed by ordinary linear differential equations with constant coeff- 
cients. They represent, in an idealized fashion, such diverse physical 
entities as motors, ships, airplanes, missiles, chemical reactors, electrical 
networks, etc. 

Linear Discrete Dynamic Elements. These are governed by linear 
difference equations with constant coefficients. Such elements are pecu- 
liar to sampling systems. They are an idealization of linear numerical 
computations performed by digital computers incorporated in the 
system. (Networks of delay lines used, for instance, to process radar 
data can be interpreted similarly.) A digital computer cannot operate 
directly on a continuous signal, but only on sequences of numbers, that 
is, on a discrete signal. Thus discrete dynamic elements inherently 
involve the operation of sampling, namely the conversion of continuous 
signals into discrete signals. 


Sample-and-Hold Elements 

These represent, again in an idealized fashion, the operation of 
“sampling’”’ various types of information (measuring a radar return 
pulse, sampling the values of any continuous signal, etc.) and then ap- 
plying the sampled information through a ‘‘smoothing”’ filter to the 
input of some continuous dynamic element. A sample-and-hold ele- 
ment is necessary whenever the output of a discrete dynamic element is 
to be applied to the input of a continuous dynamic element, since the 
output of a discrete dynamic element can only be observed through the 
process of sampling it (that is, instructing the digital computer to 
“write”? out numbers at specific instants of time) and the input of a 
continuous dynamic element cannot be a discrete signal. Sometimes a 
sample-and-hold element is inserted in the system because, paradoxi- 
cally, sampling may improve the dynamic performance of the system 


(1, 9). 


a 
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The simplest ‘‘smoothing”’ operation consists simply of holding the 
last sampled value until the next sample is available. If the input to 
the sample-and-hold element is u(t), its output v(t) is given by 


vi +r) Oc (2) 


where t, (k = 0,1, 2, ---) is the k*® sampling instant and r is the time 
elapsed since the last sample was taken. 

The manner in which the sampling instants are to be selected will be 
discussed in detail in the next section. 

Imperfect Hold. \n some cases, it is difficult to maintain the output 
of the sample-and-hold element exactly at the sampled value as required 
by Eq. 2. Instead of (2), the input-output relations of the sample-and- 
hold element may take the form 


v(t, + = (exp r/7,)u(ti), O< —h (2-A) 


where 7, is the time-constant of the sample-and-hold element. Of 
course 7, should be as large as possible compared to the average time 
interval between successive samples. 

The modification represented by (2-A) can be easily incorporated in 
the general analysis which follows and therefore need not be discussed 
further. 

Noninstantaneous Sampling. In this case, the sample-and-hold ele- 
ment transmits the input signal without change for U seconds following 
t.; thereafter, the output is held constant at the last sampled value of 
the input. Therefore in this case we have 


u(t, + T) 


u(t), 0<7rsU 
u(t, +- — ty. 


Better Smoothing. Sometimes the smoothing provided by the 
sample-and-hold element is inadequate. In such cases, additional 
smoothing may be provided by cascading a linear dynamic element 
(say, m cascaded integrators) after the sample-and-hold element. This 
modification obviously requires no special consideration. 


Interconnecting Elements 

These represent idealized operations on the outputs of the dynamic 
elements, such as multiplication by a constant, differentiation, as well 
as linear combinations of these operations. Differentiation can occur 
only in connection with the output of continuous dynamic elements. 
See in this connection, Remark (a, ii), p. 413. 

The preceding classification of elements is valid also in the case of 
nonlinear sampling systems. The latter are characterized by the fact 
that the interconnecting elements have nonlinear input-output relations ; 
in other words, the principle of superposition is not valid. 
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To illustrate the manner in which the above elements are inter- 
connected, consider the hypothetical feedback system shown in Fig. 1. 
The system consists of two continuous dynamic elements (CDE), one 
discrete dynamic element (DDE) and two sample-and-hold elements 
(SHE). The boxes with the symbols K, d/dt, and = denote intercon- 
necting elements corresponding to multiplication by a constant, differen- 
tiation, and summation. The rest of the interconnections are indicated 
simply by the lines connecting the various boxes. The quantities x; 
(the output of the discrete dynamic element) and xs, x» are discrete 
signals. All other signals are continuous. The system input is denoted 
by r(t) and the system output by x,(t). 


r(t) 


Fic. 1. A typical sampling system. 
DESCRIPTION OF THE SAMPLING OPERATION 


As mentioned in the preceding section, sampling operations occur in 
discrete dynamic elements as well as in sample-and-hold elements. 
Let us consider a single sampling operation. The kt* sampling instant 
is denoted by & (k = 0,1,---). The interval between successive 
samples 

(3) 


is called the k” sampling period. The various sampling operations may 
be characterized as follows: 


Conventional Sampling. 17, = const. for all k. 

Nonsynchronous Sampling. Vf ty, t.’ denote the instants at which 
two different sampling operations occur, then ¢,’ = & + U for all , 
where U is some positive constant. 

Multiple-Order Sampling. The sampling period is a periodic func- 
tion of k; in other words, 7, = 7.4, where q is a positive integer. 

Multi-Rate Sampling. If T,, T;’ denote the sampling periods of 


CDE, | DDE | SHE, 
Xi x, 
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two different sampling operations, 7), = const. and 7°,’ = const. for all 
k, but Ty’. 

Noninstantaneous Sampling. This can occur only in connection 
with the sampling operations of the sample-and-hold elements, whose 
input-output relations in this case are given by (2-B). 

Random Sampling. Both t, and 7, are random variables. 

The various types of sampling operations are depicted schematically 
in Fig. 2. 


kT 


2 
t' =kT'#U (T#T") 
(a) CONVENTIONAL 


SAMPLING 
(d) MULTI-RATE 
SAMPLING 


ts (e) 
SAMPLING 


SAMPLING 


they tk + Ty 
(Ty = RANDOM VARIABLE) 


lh T T, To T; 
tots fo te ty ty 


(c) MULTIPLE-ORDER 
SAMPLING (f) RANDOM SAMPLING 


Fic. 2. The sampling operation. 


Notation Convention. The following convention is used to denote 
the change in various quantities as a result of the sampling operations: 


If x(t) is some quantity whose value changes due to a sampling 
operation, and if ¢, is the sampling instant, then we denote by x(tx) 
the value of x(t) just before and by x(t,+) the value of x(t) just after the 
sampling operation. In other words, the sampling operation is thought 
of as taking place in an interval t, < t < t + ¢, where eis an arbitrarily 
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small positive number. Note that this convention is in strict ac- 
cordance with the inequalities used in (2), (2-A), and (2-B). 


MATHEMATICAL DESCRIPTION OF DYNAMIC ELEMENTS 


The mathematical description and analytic study of a dynamic 
system is greatly facilitated by focusing attention on the concept of state. 


Intuitive Definition: The state of a dynamic element is a set of 
numbers (called state variables) which contain as much informa- 
tion regarding the past history of the element as is required for 
the calculation of the entire future behavior of the element. 


The evolution of a dynamic system through time may be visualized 
as a succession of state transitions. ‘When the system is linear, these 
transitions are linear transformations of the state. The restriction to 
linearity greatly simplifies the analysis. 

The state variables of a sampling system fall into three categories: 

The state variables of the continuous dynamic elements. 
The state variables of the discrete dynamic elements. 
The state variables of the sample-and-hold elements. 


the last two cases, the state transitions may be thought of as 
taking place at the sampling instants only; in the first case, the state 
transitions occur continuously. 

We now derive the equations governing state transitions. These 
derivations are quite elementary; however, they require a certain 
amount of careful “bookkeeping.”’ In the next section, we show how 
the transition equations of the entire system are obtained by combining 
the transition equations of the various dynamic elements. The manner 
in which the equations are combined is determined by (a) the inter- 
connections of the various elements in the system, (+) the various 
sampling operations in the system. This completes the analysis; any 
information regarding the behavior of the system can be obtained from 
the transition equations. 


State of Continuous Dynamic Elements 


Since these elements are assumed to be described by ordinary differ- 
ential equations (linearity is immaterial for the moment), they may be 
simulated by means of an analog computer. ‘This simulation, on paper 
or in reality, may be done in a variety of ways; in each case, exactly 
integrators will be required, y being the sum of the orders of the differen- 
tial equations governing the various continuous dynamic elements. 
The outputs of the y integrators at time ¢t are denoted by the y-vector 
(y X 1 matrix), 

x1 (t) 


x(t) = | |, (4) 
x(t) 


af 
| 
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In the absence of inputs, the signals generated in the analog com- 
puter after t = t) depend only on the voltages stored on the integrators 
at ¢ = ty (that is, on x¢(to)) ; thus x°(¢) is clearly a correct definition of the 
state of the continuous dynamic elements (13). 

In Fig. 1, the state variables of the first continuous dynamic element 
are denoted by x1, ---, x4; those of the second by x5, x6. Thus y = 6. 

Remark a. It can be easily shown from the elementary theory of 
linear differential equations that 


(i) Any output of a continuous dynamic element is a linear function 
of the state variables. In particular, the output may always be chosen 
as one of the state variables, which is what was done in Fig. 1. 

(ii) Only those derivatives of outputs of continuous dynamic ele- 
ments are ‘“‘admissible’’ which can be expressed as a linear combination 
of the state variables of and the inputs to the dynamic element. Thus 
if v (output) and uw (input) are related by 


dy d?—y d%u 


dtp + +--+ + bu =a, 


then it follows easily that dv/dt, ---, d?~*v/dt»-« are admissible deriva- 
tives of v. 


State of Discrete Dynamic Elements 


As mentioned, these elements can be thought of as idealizations of 
numerical computations performed by digital computers incorporated 
in the system. Excluding programming instructions, fixed parameters, 
etc., suppose that 6 numbers must be stored in the computer at time ¢ 
in order to be able to perform a complete cycle of computations. These 
6 numbers are evidently the state variables of the discrete dynamic 
elements. They are denoted by the 6-vector (6 & 1 matrix) 


Xy4r1(t) | 
= | |, (5) 


In Fig. 1, the state variables of the discrete dynamic element are 
denoted by x7, xs, x» Thus 6=3. 


State of Sample-and-Hold Elements 


From (2) we see that the state of a sample-and-hold element is 
identical with its output. Assuming there are ¢ sample-and-hold ele- 
ments in a given system, we denote these state variables by the o-vector 
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i 
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(o X 1 matrix) 


Xy4041(t) 
= = | (t) 
| 


| 

which represents the state of the entire sampling system. 

In Fig. 1, the state variables of the sample-and-hold elements are 
X10, X11. Thus = 
State of Entire System 

By combining the state variables defined by (4), (5), and (6), we 
obtain the n-vector (n X 1 matrix) 
| 
x,(t) x'(t) J 


which represents the state of the entire sampling system. The integer 
= y¥ + 64+ a is called the order of the sampling system. 


x(t) = 


Transition Equations of Sample-and-Hold Elements 
The state transitions take place (discontinuously) only at sampling 
points (see (2)): 


(ty) + (te), (4 + 6 + 1, n) (7-1) 


where the s,; are constants; ¢, refers to the sampling instant of the 7 
sample-and-hold element (which may or may not coincide with the 
sampling instants of the other elements). Equations 7 express the fact 
that, in general, the input to a sample-and-hold element is a linear 
combination of the state variables of the continuous and discrete dy- 
namic elements as well as the system input r(¢). 


Transition Equations of Discrete Dynamic Elements 
Here 


Ville tr) = dijxj(te) + di ngir(te), 


j=l 


where the d;; are constants and 7; is the time required for the computa- 
tions indicated by (8) to be completed by a digital computer; ¢, is the 
instant of time at which the variable quantities entering on the right- 
hand side of (8) have been sampled (¢, may be the same for several of 
the Eqs. 8-1). 
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If the results of the computations indicated by (8-i) are not to be 
called for (as must be the case in a well-designed system) until more 
than 7; seconds after the sampling instant, ¢,, then we may just as well 
assume that the state transitions take place instantaneously at ¢,, in 
other words, the left-hand side of (8-7) may be denoted by x,(t,+) (see 
“Notation Convention,” p. 411). Often the computing time 1; is very 
short compared to the system dynamics and may be disregarded 
altogether. 


Transition Equations of Continuous Dynamic Elements 
Let us define, in conjunction with the analog computer setup intro- 
duced in a previous section (see p. 412), the constants 


a;, = coefficients with which the output of the j integrator is fed 
back to the input of the zt" integrator. (See Fig. 3.) 


Fic. 3. An analog computer simulation. 


In the absence of any external inputs to the continuous dynamic ele- 

ments (13) we have by inspection of Fig. 3, 

j=l 

This equation is a quantitative description of the state transition 
x,(t) > x,(t) + dx,(t) due to the infinitesimal time change ¢t — ¢t + dt. 
To calculate the state transition for a finite change in time, it is neces- 
sary to solve the system of differential equations (9). We first cast (9) 


dx dx 
nat n Yat 
@) 
: 
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into the standard vector-matrix notation: 


dx‘(t) 
Ax*(t) (10) 


where x°(t) is a vector defined by (4) and A isa (y X y) constant matrix 
whose elements are defined above. 

It is well known (14) that solutions of (10) starting at any time ¢) and 
observed at any time ¢ are given by 


x°(t) = (exp (t — to)A)x*(to) = — to)x*(to) (11) 


where 


exptA = t*A*/k! (A® = unit matrix). (12) 
k=0 

The infinite series (12) converges for all finite values of tf. **(¢) is 
called the transition matrix (11) of the continuous dynamic element. 

The form of solution given by (12) is quite convenient if ¢ is some 
fixed constant (as will frequently be the case in the examples) because 
it does not require finding the roots of the characteristic equation of the 
matrix A. On the other hand, it may be inefficient to use (12) for 
finding &°(t) for a large number of values of ¢. 

A convenient way of obtaining a closed expression for 6°(t) is the 
following. Let 


gi;(s) = transfer function from input of j‘* integrator to output of 
z** integrator in the analog-computer setup of Fig. 3. 
Then 
gu(S) 
= | = £[G(s)] (13) 
£y1(S) £17(S) 
where the inverse Laplace transformation £~! is to be applied separately 
to each element of the matrix. The calculation of transfer functions 
and their inverse Laplace transforms can be carried out by well known 
methods of linear system engineering. In the examples which follow, 
calculation of &¢¢(t) by means of (13) will therefore not require additional 
explanation. 
The following important properties of &(t) should be noted: 
+ to) = (t,)b (te) for all ty, te. (14) 
This follows directly from (12). Therefore, in particular, we have 
(0) = unit matrix; ((t))-! = 2). (15) 


To complete the derivation of the transition equations for the con- 
tinuous dynamic elements, consider now the inputs to these elements. 
As mentioned previously, the external inputs to these elements are 
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applied through sample-and-hold elements. (There may be inputs 
from one continuous dynamic element to another but these are im- 
plicitly taken care of in the definition of the A matrix in (9).) 

When the sampling operations are nonsynchronous, multiple-order, 
multi-rate or random, the output of the sample-and-hold element is a 
constant between successive sampling instants. Therefore the inputs 
to the continuous dynamic elements will be a series of steps. 

Let us replace, for a moment, the j*® sample-and-hold element by 
means of a generator whose output at time ¢ = ¢& is a unit impulse. 
The effect of the impulse is to change instantaneously the state of the 
continuous dynamic elements from x°(t,) to x°(t.+) = + Ax*(t,). 
Let us denote the change by the y-vector (y X 1 matrix) 


fi = Ax*(t,). 
In view of the linearity, the effect of the unit impulse on the succeeding 
time-variation of the state variables is expressed by 
x°(t) = (exp (t — &)A)f, t >t. 
If instead of applying a unit impulse, the generator in the place of the 7 
sample-and-hold element applies a step of magnitude x;(t,+), we get, 
again using linearity, 


Ax‘(t) = (exp u)A)fidu t>t 


tk 


x; (t,*) (exp uA)fidu x;(t,+)h(t), t> tk. (16) 
0 


Using (12) and (13), the following explicit formulas are obtained for the 
vector h’(t) defined by (16): 


hi(t) = /(k + 1) (17) 
= | (18) 


We can now write down the complete transition equations for con- 
tinuous dynamic elements. Let ¢, be the time at which any one of the 
sample-and-hold elements changes it state. Let ¢,4; be the instant at 
which the next change of state of any one of the sample-and-hold ele- 
ment occurs. Then 


+ 7) = + 


(@=1,2,---,y) (19-2) 


It remains to consider sample-and-hold elements in which the 
sampling operation is of the noninstantaneous type described on p. 406. 


\ 
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The output of such elements after t, + U (see (2-B)) is constant, so 
that the previous analysis applies. During the interval t, <¢t =t, + U, 
the output of the sample-and-hold element is equal to the input. The 
input to the sample-and-hold element must consist of combinations of 
state variables of continuous dynamic elements and the system input 
r(t); no state variables of discrete dynamic elements can appear. 

Let the x,(t) be the state variable of the only sample-and-hold ele- 
ment of the noninstantaneous type described on p. 409. By (7), and 
bearing in mind the above remarks, we have 


Xn(t) = SurXi(t) + +--+ + +  (7-n) 


during the time interval 4, <¢ = % + U. We now form a matrix A 
by defining its elements as 
ai; = ai; + 


Similarly, we define a vector f = s,,,4:f". Then the transition equa- 
tions in the interval t, <¢ =t4, + U are obtained by solving the 
differential equation 


dx*/dt = Ax: + fr(t) (20) 


which yields 


x°(t, + 7) = (exp tA)xe(t,) + r(u)(exp(r — u)A)fdu, 
0 


An analogous procedure may be employed when there is more than 
one noninstantaneous sampling element. The tedious, but elementary, 
details of setting up equations in those cases are left to the interested 
reader. 

SYSTEM TRANSITION EQUATIONS 


The system transition equations can now be written down very 
simply. The state transitions of each one of the three types of elements 
are represented by a matrix acting on the combined state variables. By 
multiplying together the various transition matrices, we obtain the 
transition matrix for the system for any two successive values of time. 

We assume first that r(t) = 0. 


Sample-and-Hold Elements 
Let S; (¢ = y + 6+ 1, ---, m) be the matrix obtained by replacing 
the row of the X m unit matrix by Si2, «++, Si, 748,90, 0. 


(See (7-1).) 
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| 

| 

| 

1 
| 
| 
| 


1 
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If the state x(t,) of the over-all system is multiplied by the n x n 
matrix S;, then all states remain the same except the z** which changes 
in accordance with (7-7). Note also the relations 


SiS; = SiS; and = 


Discrete Dynamic Elements 

Let Dj;...(4, Jj, =y+1,---,y +4) be the matrix obtained by re- 
placing the row of the m X m unit matrix by di, ---, di 445, 0, 0; 
replacing the row by dj, (See (8-2), 8-7)): 


row (23-ij- ++) 


— j* row 


If the state x(t) of the over-all system is multiplied by the n x n 
matrix D,;... then all the states remain the same, except x,(t), 
x(t), which undergo a transition in accordance with (8-7), (8-7), 
all at the same time t;+ = t;+ = --- 


Continuous Dynamic Elements 


Let @(7) be the matrix defined by: 
( | | 
hyvt8+1 (7) 
e«(r), O , : : first y rows 


(24) 


a 

I 
419 

0 

0 1 0 

S; = (22-1) 

1 0 O 

| 

Djj... = | 
dj, +++ djy | Gj | 

| 

|--- 

0 

| | 
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If the state x(t) of the over-all system is multiplied by the n X n 
matrix (7), then the states of the discrete dynamic elements and 
sample-and-hold elements remain the same and the states of the con- 
tinuous dynamic elements undergo a continuous transition as a function 
of +. The argument of #(7) may have any value in the interval 
0 <7 S bey: — tk, where ¢, is the last sampling instant and ¢,4; is the 
next sampling instant. 

The system transition matrix, for any interval of time, is formed by 
multiplying together the individual matrices S,, Dj;... and ®(7) as re- 
quired by the various sampling operations taking place in the system. 


Effect of System Input 

We now remove the restriction r(t) = 0. For simplicity, we as- 
sume that r(t) can affect the system only through a sample-and-hold 
element. This will exclude minor complications due to the convolution 
integral in (21). In accordance with the previous generalized notation, 
we now define the n-vectors 


which correspond to the coefficients of r(¢) on the right-hand sides of 
(7-2) and (8-2). These vectors show the effect of r(t) on the state 


transitions. 


Example 1. System Defined by Its Transition Equations 
The over-all transition equations of a given system have the form: 


= g)SioP(T 4) 11 (te) 
+ 4) (dy + Sii)r (ty). (26) 


The question is: What kind of a system is the one whose transition 
equation is given by (26)? The answer can be obtained by systemati- 
cally examining the various terms in (26). 


(a) There are eleven state variables. The first six state variables 
belong to continuous dynamic elements, the next three to discrete 
dynamic elements, the last two to sample-and-hold elements. 


0 | 0 
0 0 
d,; = s; = 0 (25) 

0 

0 0 

= 


May, 1959.] THEORY OF SAMPLING SYSTEMS 421 


(6) Each sampling operation has the same period 7 = 74 + Ts; 
the sampling operations in connection with the discrete dynamic element 
and the second sample-and-hold element (whose output is x;;) occur si- 
multaneously at ¢,; the sampling operation in connection with the first 
sample-and-hold element (whose output is x10) occurs at t’ = th + Ta. 
Thus the sampling is of the nonsynchronous, fixed-period type. 

(c) All information regarding the input r(t) is obtained at t,. The 
first sample-and-hold element samples only the state variables of the 
continuous and discrete dynamic elements. This may happen, for 
instance, in order to utilize the computations performed by the dis- 
crete dynamic element (which may be assumed to require less time than 
T 4 for completion) before the next sampling point ¢,41. 


It is easy to see that (26) may represent the operation of a system 
such as Fig. 1. 


f(kT) m(t) 


e(t) 


Fic. 4. Example 2. 


SOME ILLUSTRATIVE EXAMPLES 


To show how the formalism of the preceding sections is applied 
to the derivation of transition equations of systems containing various 
types of sampling operations, one example for each type of sampling 
will be given in this section. 

In each example, the derivation of the transition equations requires 
essentially four steps: 


Step (i). Selection of the state variables. 

Step (ii). Calculation of the various matrices and vectors defined in 
the two previous sections. 

Step (11). Derivation of the transition equations for the various 
sampling operations and sample-free time-intervals. 

Step (tv). Combination and simplification of the transition equations. 


Example 2. Conventional Sampling System 


A typical system of this class is shown in Fig. 4. The system con- 


i 

é 

: 

q 

: 

7 

Be 

: 


422 R. E. Karman anp J. E. Bertram (J. F. 


tains a discrete dynamic element (for improving the stability of the 
system). ‘The sampling operations in the discrete dynamic element 
and the sample-and-hold element are performed synchronously. (This 
is not always necessary or even desirable in practice. See Kranc (6).) 
Thus sampling occurs at the instants: 


Discrete Dynamic Element: 0, T, RT, 
Sample-and-Hold Element: 0, T, kT, 
(i) The continuous dynamic element in Fig. 4 is described by the 
differential equation 
(27) 


e(t) r(t) 


which relates the output c(t) of the system to the control signal m(t). 
The discrete dynamic element is described by the difference equation 


f(RT) = ace(kT) + ae((k — 1)T) — bi f((R —1)T). (28) 


The input signal e(t) to the discrete dynamic element is continuous. 
This signal is sampled and only its samples affect the state transitions 
in the discrete dynamic element. The output f(k7) of the discrete 
dynamic element is a discrete signal, that is, it is available only at 
time kT. 

The state variables for the continuous dynamic element are selected 


as follows: 
xi(t) _ c(t) 


CDE | 
x, Xo (t) x(t) | 
| | 
| 
| 
DDE 

| 

| 

| | 

| 

Fic. 5. Example 2 (details). 
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With this choice of the state variables, the block diagram for the 
continuous dynamic element governed by (27) is shown in Fig. 5. This 
block diagram shows also how Eqs. 27 could be simulated by means of 
an analog computer. 


To define the state variables for the discrete dynamic element, (28) 
is written in the form 

F(RT) = ace(RT) + (a1 — aobi)y((Rk — 1)T) (30) 
WRT) = e(RT) — biy((k — 1)T). (31) 
By eliminating y(k7T) and y((k — 1)7) from (30) and (31), it follows 
easily that Eqs. 30 and 31 are equivalent to (28). From this, it is clear 
that the state of the discrete dynamic element at time kT is given by 

y((k — 1)T). Hence let 
x"(kT) = [xa(kT)] = Ly((k — 1)T)]. (32) 
Finally, the state variable of the sample-and-hold element is evidently 


given by: 
x"(kT) = = [m(kT)] = [f(RT)]. (33) 


(it) The transition matrix for the continuous dynamic element is 
calculated according to the discussion given on p. 419; the end result is, 
using (24), 

r—1i+e"’ 
| 


0 


1 


| 
| 
| 
| 
| 


Using the fact that e(t) = r(t) — x,(¢) and the notation of (32), we can 
write (30) and (31) in the form 


x3((R +1)T) = — xi(RT) — bix3(RT) + r(kRT) 
f(RT) = (ay — aoby)x3(RT) + (RT) — x,(RT) ]. 


Hence, from (23) and (25) 


= 

0 0 1 | 

a 

1 0 0 0 

| | 

0 0 0 


If Eq. 30 is changed similarly, using the notation of (33), we get 


x(RT*+) = — aoxi(RT) + (a, — aobi)x3(RT) + aor(RT). 


Hence, from (22) and (25), 


(iii-iv) To obtain the transition equations for the over-all system, 
we must examine the sequence of the various transitions. Consider 
the time ¢ = kT. First, the new value of the control signal x,(kT) is 
computed and we have 


x(kRT+) = Syx(kT) + sar(kT). 


Next, the data for the computation of the state transition of the dis- 
crete dynamic element are collected and the new state x; is calculated. 


x(kT++) = Dyx(kT*) + der (RT) 
= DS x(kT) + (ds + 8a)r(RT). 


Actually, the new value of x; may not be available at kT++ but only 


(39) 
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somewhat later; this is immaterial, however, as far as the transition 
equations are concerned. 

During the interval RT < ¢ = (k + 1)T the only transitions taking 
place are those of the state variables of the continuous dynamic element. 
Hence, 
x(kT + 7) = + (ds + 8a)r(RT)], O< = T (40) 

= W(r)x(kT) + v(r)r(k7), O<raT (41) 


where 


V(r) =6(7)D3S,= 


1 —ao(r (r —1+e7") (a1 


> 


-—i+e") 


v(r) = &(r)(ds + 84) = 


Equation 41 describes the behavior of the system at any instant of time. 
Remark b. By letting r = T in (41), we have a stationary (constant- 
coefficient) difference equation 


x((k + 1)T) = W(T)x(kT) + v(T)r(kT) (43) 


relating the state variables at succeeding sampling instants. This is 
equivalent to the result usually obtained by means of the z-transform 
method. 

Remark c. The last column of the matrix (7) contains only zeros. 
This shows that the value of x, need not be known at the k* sampling 
instant in order to find the state vector at the (k + 1)** sampling in- 
stant. Thus x, may be dropped in Eqs. 41 or 43. If x4 had been 
dropped before obtaining the final result, however, one could not have 
expressed (7) as a simple product of matrices. 


0 
| 
| 
—1 0 | —b, | 0 
| 
0 | ay, — ab; 0 
: 
and 
ao(1—e-") 
|- - - - - - - (42-B) 
~-------- 
ae 

a 

-- 
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Remark d. \n deriving (41) and (43), no real use was made of the 
fact that 7, = T = const. Hence these equations are valid also when 
the sampling period is not a constant. In that case the transitions 
between successive sampling instants are governed by the nonstationary 
(varying-coefficient) difference equation 


= + (te) (44) 


where ¢, denotes the kt" sampling instant and T, = ti41 — tk. 


Example 3. Multiple-Order Sampling 


We consider the same system as in Example 2. We assume again 
that the sampling operations are synchronous, but now they take place 
at the instants 


=0, To+T, Tot: +7,-1.=T, kT, kT+T,, 
In other words, the sampling period 7, is a periodic function of k with 
period 
(i-iii) This example is evidently a special case of the nonstationary 
system (44) derived in connection with the previous example. 
(iv) Because T, is periodic, we may relate x(kT) to x((k + 1)T) 
by means of the stationary difference equation : 
x((k + 1)T) = W(To, Ty-1)x(RT) 
+ vi(To, To-1)r (RT) 
+ T,-1)r (kT + To) + 
+ + To + + Ty-2) (45) 


where 
W(To, Te-1) = V(T 9-1) V(T¢-2) Y(To) 


vil To, 1) )¥(Ty-2) v(T») (46) 
vq( 74-1) v(7,-1), 
where V(r), v(7) are given by (42-A-B). 


Example 4. Nonsynchronous Sampling 
The system is shown in Fig. 6. The sampling operations occur at 
the following instants of time: 
Sample-and-Hold Element 1: 0, kT, --- 
Sample-and-Hold Element 2: U, T+ + U,::: 


= 
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For simplicity, it is assumed that 7 > U > 0. 


(i) The definition of the state vector is apparent from Fig. 6. 
(ii) The transition matrices for the two sample-and-hold elements 
are given by 


Fic. 6. Example 4. 


while the transition matrix of the continuous element is 


0 


| 

| 
— | - - 

| 0 

| 


427 4 
q 

| | | 

1 0 | 0 1 0; 0 | 0 

0 1; 0; 0 0 1! 01 0  . 
S,= |-----1---1---| and S,= |----!---!--- 
x x 

3 

CDE, 

4 

| | 

CDE, 

| | 

| 

| 

Xo 

| 

£44 

0 


428 R. E. KALMAN Anp J. E. Bertram 


From Fig. 6 we have also 


(iii) At time t = k7, there is a discontinuous change in the value of 
x3 which is expressed by 


x(k7+) = S3x(kT) + sar(RT). 


(47) 


During the time-interval RT < t = kT + U no sampling operations 
occur and the state transitions are continuous. Thus 


x(kT + 7) = (r)x(kT+), O< 7 =U. (48) 


At time t = kT + U there is another discontinuous state transition 
due to a change in the value of x,: 


x(kT + Ut) = S,xx(kT + UV). 


During the time-interval RT + U <t = (k + 1)T, the state transi- 
tions are again continuous: 


x(kT + U+ = + Ut), O< r ST-U. (50) 


(49) 


(iv) Equations 47-50 describe the behavior of the system at every 
instant of the time interval RT <¢ = (k + 1)T. Confining attention 
to the transition from x(k7) to x((k + 1)T), we get the following sta- 
tionary difference equation by combining (47-50) 


x((k + 1)T) = W(T, U)x(kT) + v(T)r(kT) (51) 


where 


W(T, U) = &(T — U)S@(U)S;3 = 


1 


U. F. 1. 
0 
0 
8; = |--- and 8, = 0. 
1 
| 0 0 
| 
0 - 1 | 0 0 
1 - U 0 0 
| 
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and 


U(1 e~(T-U)) 
v(7, U) = (T — U)S,®(U)s; = |--------- 


Remark e. Notice that in (51) the state variable x; (but not x,4) is 
superfluous. 

Remark f. As in Example 3, the derivation of the above equation 
does not make use of the assumption that 7, U are constants. Letting 
these quantities be a function of k we get the nonstationary difference 
equation 


X(te41) = (ty) + v(7%, U;,)r (te) (53) 


where ¢, denotes the k** sampling instant of SH7E,, and 


— ty = > > 0. 


Example 5. Multi-Rate Sampling 

Consider again the system shown in Fig. 6. In this example, the 
sampling operations have periods 74 = 7/3 and T, = T/2, respec- 
tively. Thus sampling occurs at the instants 


Sample-and-Hold Element 1: 


kT, 7/3, kT + 27/3, (Rk +1)T, - 


Sample-and-Hold Element 2: 


kT, RT +T/2, (k+1)T,-* 


(i-ii) The transition matrices of the two sampling operations and of 
the continuous dynamic element are the same as in Example 4. 

(iii) At time t = kT, both sample-and-hold elements have a dis- 
continuous state transition : 


x(RT+) = + sar(RT). 


During the interval kT < t = kT + T/3 the state transitions are 
continuous 


x(kR7T + r) = ®(r)x(RT+), O< = T/3. 
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At time k = kT + 7/3, 
x(kT + T/3+) = Sx(kT + T/3) + sar(RT + T/3). 
During the interval RT + 7/3 <t SkT + T/3, 
x(kT + T/3 + 7) = O(r)x(RT + T/3+), O< 7 ST/6. 
At time ¢ = kT + T/2, 
x(kT + T/2+) = Sx(kT + T/2). 
During the interval kT + 7/2 <7 SkT + 27/3, 
x(kT + 7/2 + r) = O(r)x(RT + T/2+), O< S=T/6. 
At time t = kT + 2773, 
x(kT + 27/3+) = S3x(kT + 27/3) + sar (RT + 27/3). 
During the interval RT + 27/3 <t S (k + 1)T, 
x(kT + 27/3 + 7) = O(r)x(RT + 2T/3+), O< 7 ST/3. 


(iv) Since the pattern of sampling operations repeats with a com- 
mon period 7, by combining the preceding equations we can obtain a 


stationary difference equation relating x((k + 1)T) to x(kT): 


x((k + 1)T) = Wx(kT) + vir(RT) + vr(kT + T/3) 
+ var(kT + 27/3) (54) 


where 
W = 
= 
ve = 
v; = &(T/3)s;. 


Remark g. It should be carefully noted that the derivation of the 
stationary difference Eq. 54 was possible only because the ratio of 
sampling periods was rational. In fact, if 74/T,» = p/g (where p, g are 
relatively prime integers), then the sampling patterns have a least 
common period g7’4 = pT, = T and therefore it is possible to relate the 
state x(k7) to x((k + 1)T) by means of a stationary difference equa- 
tion. When 74/7, is irrational or, what is practically the same thing, 
the integers ~, g are large, then a very large number of intervals have to 
be considered in order to write down the over-all stationary transition 
Eq. 54. 


(55) 
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Remark h. Alternately, it is possible to write down a nonstationary 
difference equation of the type 


x((k + 1)T) = V(k)x(kT) + v(k; r(t)) (56) 


where T is not the least common period of 7, and T,. In this case, the 
matrix V(k) and the vector v(k;7(t)) must be calculated separately for 
each interval RT < t = (k + 1)T by considering the various sampling 
instants and the various sampling-free subintervals as was done above. 
When 74/7 z is irrational, that is, when the least common period of 7°, 
and 7’, is infinite, the matrix ¥(k) will have coefficients which vary in 
a quasi-random (nonperiodic) fashion with k. Therefore under such 
circumstances one may expect the sampling system to exhibit a be- 
havior which would be intuitively interpreted as ‘‘random.”’ 


Example 6. Noninstantaneous Sampling 

This system is shown in Fig. 7. The sampling operation is as de- 
fined by (2—B.) 

(i) The definition of the state variables is apparent from Fig. 7. 


x x 
CDE 


Fic. 7. Example 6. 


element. During the time interval RT < ¢t = kT + U, the sample-and- 
hold element transmits all signals without modification. Thus 


x2(t) = r(t) — x(t). (57) 


During this time interval, the continuous dynamic element is governed 
by the differential equation 


dx,/dt = Gx, + r(t). (58) 


Integrating (58) and combining it with (57), we obtain the transition 
equation 


| 

| 

r(t) o—+(S)—+ SHE | 

(ii-i11) There are two transition matrices for the continuous dynamic [ae 
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x(kT + 1) = + fe — + + 1), 
0<7rsU (59) 


where 


(60) 


= i--- = |---|. 
0 | 1 


During the time interval kT + U =t = (k + 1)T, the output of 
the sample-and-hold element remains constant and equal to x2(kT + U); 
the continuous dynamic element is governed by the differential equation 


dx,/dt = Gx.(kT + U). 
From this we find by inspection the transition equation 


x(kT +U+7) = @(r)x(RT+ VU) OS7ST-U (61) 


where 


-|. (62) 


| 

(7) = |---1- 

| 

| 


J 


(iv) Combining these results, we obtain an expression for x((k+1)T) 
in terms of x(kT) as follows: 


x((k + 1)T) = ¥(T, U)x(kT) + v(T, U, r(t)) (63) 


where 


and 
v(T, U, r(t)) = &(T — VU) [sarcer + U)+ 


It should be observed that (63) is not a true difference equation (but 
rather a difference-integral equation) because of the complicated fashion 
in which the system input r(¢) affects the state transitions. 


| 
ear 0 
= |---1--- 
— er 0 
| 5 
| | 
ar 
1 
—U)); 0 
(64) 
— 
| 
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Example 7. Other Applications 

There are many other problems in the analysis of nonstationary 
linear dynamic systems which can be handled by our methods. We 
mention two further examples: 

(1) In discussing the solutions of Mathieu’s equation with a square- 
wave parametric excitation, that is, the equation 


d*x/dt? + adx/dt + b(1 + ¢ sgn (sin wt)) = 0 


where a, b, c, w are positive constants and sgn(0) = 0, sgn x = x/|x| if 
x #0, Pipes (15) has used a method which is conceptually identical 
with ours. 

(2) In a recent paper, Desoer (16) has discussed a dynamic be- 
havior of networks which contain periodically operated switches. 
Desoer derived the equivalent of formula (12), apparently regarding this 
as a new result. In applying our method to Desoer’s problem, very 
considerable conceptual and analytic simplifications may be obtained. 
We omit the details. 

STABILITY 


Probably the most important task of system theory is to answer the 
question, “Will a given dynamic system be stable?”’ The first sig- 
nificant result of the theory of (conventional) sampling systems was the 
derivation of a stationary difference equation governing the behavior of 
the system between successive sampling instants; the question of 
stability can be answered directly by examination of the difference 


equation. 
Following a suggestion of Zadeh (17), the concept of stability is 


formulated by the following 


Definition—A linear system is said to be stable if and only if 
every bounded system input signal produces bounded varia- 
tions in all state variables. 


This definition applies to linear systems regardless of whether they are 
governed by differential or difference equations or whether they are 
stationary or nonstationary. 

In a linear system, stability depends only on the transition matrix 
of the system; stability cannot be brought about or destroyed by a 
particular choice of the initial state or the system input signal. It is not 
difficult to show that the above definition implies. 


Theorem 1. A stationary linear system 1s stable if and only if 
every element of V™ tends to 0 as N — @ [every element of 
W(k)W(k + 1)---W(k + N) tends to 0 uniformly* in k as 


N 


‘It can easily be shown that if ¥(k)W(k + 1) Wk + N) does not approach 0 uniformly 
as N— «, then for some bounded system input signal some state variable will become un- 
bounded. 


| 
: 
| 
: 
| 
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In the above theorem, V() is the transition matrix of the system 
over some well-defined, constant interval of time. 

One can “‘experimentally” verify from Theorem 1 the stability of a 
given system by multiplying together the transition matrix a large 
number of times. This is admittedly a brute-force procedure, but, 
because of the simplicity of matrix multiplication, very well suited for 
machine computation. Moreover, in the case of a nonstationary 
system, this method may be the only way of checking stability since no 
explicit necessary and sufficient stability conditions for nonstationary 
linear systems are known at present. (There are, of course, many 
sufficient and many necessary conditions.) 

It is fortunate that in many sampling systems the transition matrices 
over sufficiently long intervals of time are constant ; see Eqs. 43, 45, 51, 
54 and 63. (For cases where the transition matrix is not constant, see 
Eqs. 44, 53, 56 and Remark g). For a constant transition matrix, the 
explicit stability conditions are as follows: 


Theorem 2. A stationary linear system is stable if all n (possibly 
complex) zeros, z;, of the polynomial det (VW — zl) satisfy 


<1 Ds (65) 


If, moreover, |z:| # 1 for any i, then condition (65) ts not only 
sufficient but also necessary. 


The proof of this theorem follows readily by reducing W to its 
Jordan canonic form. 

Condition (65) is identical with the result that the poles z; of the 
z-transform of the input-output relations of the system must lie within 
the unit circle (1). Of course, our result is much more general than the 
z-transform result because Y can be computed in many cases where the 
z-transform method is not applicable. 

Finally, it should be remarked that a stationary transition matrix 
can be obtained in those cases and in those only where the pattern of the 
sampling operations repeats in a periodic fashion. This observation 
generalizes Floquet’s theorem concerning the nonstationary differential 
equation 


dx/dt = A(t)x (66) 


where A(t) = A(t + T), that is, the system is periodic with period T. 
According to this theorem (14) the transition matrix v(t + 7; t) of (66) 
for an interval of length 7, 


x(¢+ 7) = W(t + 7; 2)x(t) 
is given by 
W(t + 7;t) = exp BT 


4 
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where Bisa constant matrix. This shows that the stability problem ina 
nonstationary but periodic linear system is essentially the same as in 
a stationary linear system provided that attention is focused on the 
transition matrix corresponding to one period. 


CONCLUSIONS 


The development of the theory and engineering application of 
sampling systems has been hampered by the fact that hitherto standard 
methods of analysis (the z-transform method and its various modifica- 
tions) are not readily applicable unless all sampling operations in the 
system are synchronous and of constant sampling period. This paper 
shows that such restrictions are not at all basic; sampling operations of 
any presently conceivable kind can be readily studied within the unified 
framework presented here. 

Sampling systems are a special class of nonstationary dynamic 
systems. As a result of the nonstationarity, whose complexity is 
directly dependent on the complexity of the sampling operations, the 
description of such systems is necessarily somewhat involved. In many 
cases the behavior of sampling systems is still stationary if attention is 
focused not on the instantaneous state transitions but on the transitions 
over a certain time-interval. The transition matrix for these intervals 
can then be used to answer the most basic question of system theory, 
namely that of stability. 

In fact, any other information concerning the system behavior can 
be derived from the transition equations. Thus, the classical results 
concerning the stability of sampling systems can be extended to cases 
where the sampling operations are not restricted to the synchronous, 
constant-sampling-period case. 
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THE FRANKLIN INSTITUTE 


MINUTES OF THE STATED MEETING 
April 15, 1959 


The Stated Meeting of The Franklin Institute was held at 8:15 P.M. in the Lecture Hall. 
Mr. Wynn Laurence LePage, President, called the meeting to order with approximately 
225 members and guests in attendance. 

The President stated the minutes of the March 18, 1959 meeting will be published in the 
April 1959 issue of the JoURNAL and will be presented for approval at the next meeting. 

The President noted that this was the occasion of the Sixth Edward G. Budd lecture, 
which lectureship was founded in 1948 in honor and in memory of our distinguished former 
member of the Institute and its Board of Managers. Until his death in November 1946, 
Mr. Budd had been a member of the Institute since 1895 and a member of its Board of Managers 
since 1927. 

The President then introduced the Budd Lecturer, Mr. Henry Belin duPont, Vice President 
of the E. I. duPont de Nemours & Company, who spoke on ‘‘The Greatest Invention of Them 
All.” 

Mr. duPont spoke of the real progress that has been made in technology, in the 135 years 
since The Franklin Institute was founded, noting Philadelphia’s rich technology history. 
He recounted the difficulties faced by early inventors, not only in getting backing for their 
ideas but in actually creating inventions, for often they could only make a part-time job of it 
while they made their living at something else. He brought out that the modern industrial 
establishment has made the pursuit of technology a full-time job, not only for a few but for 
hundreds of thousands. Moreover, the development of the modern corporation has helped 
to increase and broaden the resources. It has provided a focal point for the development of 
technology and, while the research efforts of the universities, private foundations and govern- 
ment laboratories contribute importantly to the advancement of American science and engineer- 
ing, these institutions could not exist were it not for the wealth generated by American industry. 

Following Mr. duPont's talk there was an active and most interesting question period. 
Before adjourning the meeting at 9:20 p.m., the President expressed on behalf of the Institute 
its warm thanks to Mr. duPont. 


F. JACKSON, JR. 
Secretary 
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COMMITTEE ON SCIENCE AND THE ARTS 


(Abstract of Proceedings of Stated Meeting held Wednesday, April 8, 1959.) 
HALL OF THE COMMITTEE, 
PHILADELPHIA, APRIL 8, 1959. 


Dr. G. Scumipt in the Chair. 
The following report was presented for final action: 
No. 3333: Work of John Hays Hammond, Jr., in Remote and Automatic Controls. 
This report recommended the award of an Elliott Cresson Medal to John Hays Hammond, 
Jr., of Gloucester, Massachusetts, “In consideration of his many inventions and developments 
in the fields of radio and electronics; and, particularly, for his pioneering contributions and 
leadership in the invention, design and development of remote radio control for vehicles moving 


on or under the water. 
D. S. FAHRNEY 


Secretary to Committee 


YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 135-year history, there is vital need for 
The Franklin Institute to effectively promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive, for they begin with students in 
the early grades of our elementary schools and continue throughout an individual’s 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology through education. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 
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THERMOELECTRIC ENERGY CONVERSION IN SOLIDS 
BY 
NICHOLAS FUSCHILLO! 


I. INTRODUCTION 


We are presently in the midst of a period of great change in electrical 
technology. This change has hitherto been greatest in the field of 
electronics where a host of solid state phenomena and devices are finding 
application as circuit elements undertaking diverse tasks such as ampli- 
fication, rectification, oscillation, and detection and interconversion of 
mechanical, electromagnetic and electrical impulses. To this formid- 
able list has recently been added semiconductor thermoelectricity with 
its capability for the direct conversion of heat into electricity and the 
direct electrical transport of thermal energy from one place to another. 
This development could profoundly affect modern technology which is 
of course largely based on the interconversions of heat and electrical 
power. In semiconductor thermoelectric devices, these interconver- 
sions take place without recourse to mechanical machinery, and con- 
siderable simplifications as well as other obvious advantages can be 
achieved by the elimination of dynamos, electric motors, compressors, 
and working fluids in electrical power generation, refrigeration, and 
heating. The physical effects upon which these conversions of elec- 
trical and thermal energy are based (the effects of Peltier, Seebeck, 
Thomson, and Joule) have been well known for more than a century. 
Nevertheless, it was not until the rapid development of semiconductor 
physics in the last two decades that the conditions for maximization of 
the effects concerned were properly understood. This understanding 
was achieved through basic research in semiconductor physics and has 
led in the past few years to the emergence of these phenomena from the 
physicist’s laboratory into the realm of highly promising engineering 
developments. 

After quite a short period of intensive research, considerable progress 
has been made in the development of suitable materials for thermo- 
electric engines. Already existing materials are capable of power 
generating efficiencies approaching 10 per cent over useful temperature 
ranges, and efficiencies of 30 per cent may be achievable in a few years. 
In thermoelectric refrigeration materials, progress has also been en- 
couraging. Figures of merit of 3.0 « 10-* deg-! at 300° K and tem- 
perature differences of 75° from 300°K have been obtained with single 
thermoelements. Moreover, the efficiency of thermoelectric heating 
can be made several times more efficient than ordinary electrical re- 
sistance heating due to the thermoelectric heat pumping action. Natur- 
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ally, important problems still await solution. The figures of merit and 
efficiencies of heat pumping materials are too low by a factor of two or 
three for wide-spread conversion to thermoelectric methods in refrigera- 
tion and air conditioning. For thermoelectric generation, figures of 
merit of 3.0 & 10-* deg~' cannot as yet be maintained over a sufficiently 
large temperature to enable the attainment of energy conversion 
efficiencies competitive with the heat engine and dynamo combination. 
A number of additional problems also exist. Prominent among these 
is the contact resistance and bonding difficulties, which limit 
the miniaturization of thermoelectric devices, and the economical 
arrangement of effective heat exchange with the surroundings. We 
may expect, however, considerable progress to be made in the next few 
years in solving these problems in spite of their difficulty. 


2. ENGINEERING RESULTS OF THE ZERO-ORDER PHENOMENOLOGICAL THEORY 


The basic thermoelement (also known as a thermocouple) used as a 
thermoelectric engine is shown in the accompanying figure. Here, 
n-type and p-type semiconductor legs (unshaded) are connected by 
metallic current leads (shaded) so that a current, J, flows through them 


S2, Ke, 2 


$i, Ki, a 
= 
r2, Re, A2,4 


p-type 


ri, Wr, Ar, 2 


4-44 


V=1I R 
m=R/( 


in series. Referring to either leg 1 or 2 of the thermocouple we intro- 
duce the following notation; where the first four quantities represent 
average values over the appropriate temperature range. 


Seebeck coefficient, volts/° K. 
specific electrical conductivity, ohm~! 
specific thermal conductivity, watts-cm-! deg™! 
total electrical resistance, r = 7; + 72 
total thermal conductance, k = k; + ke 
uniform cross section area of a leg 

= length of a thermocouple leg, /; = J: = 1. 


K 
k= 
A 
l 
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As n-type and p-type semiconducting materials are connected in series 
the Seebeck coefficient of the combination S = S,; — S, will be additive 
so that its magnitude |S| = |.S,| + |S.|. 7 and To are the absolute 
temperatures of the variable and controlled (thermal sink) junctions of 
the semiconductor legs of the thermoelement. This temperature differ- 
ence is caused by the reversible absorption and evolution of Peltier 
heat at T and 7». 

When used as a heat pump we may define the thermoelement effi- 
ciency by a coefficient of performance (C.O.P.) given by (see Ref. 1) 


6 =C.0.p, = Rate of absorption or evolution of heat at T _ Qr 


Total electrical power consumed Ww 


When the thermoelement is used as an electrical power generator, 
we consider the external current source replaced by an external circuit 
of load, R (shown dashed in the figure). The current, 7, developed 
from the thermoelectric Seebeck e.m.f., is fed into the external circuit 
of impedence, R, where it does useful work. The useful power delivered 
to R will, of course, depend on the impedence match ratio m= R/(r:+72). 
The technical efficiency of the generator is defined as the ratio, 


Useful electrical power, /?R, delivered 
to the external circuit 


7 Rate of heat energy absorption at the’ 
variable junction at 7° K. 


The principal results for thermocouples used as heat pumps and 
power generators with geometrically optimized areas are, as summarized 
in cases 1, 2 and 3 below provided we have: 


Ai _ N = V1 + 4Z.(T + 7»), 


and the figure of merit, 


(Si — S2)? 


VK2/02)" 


Case 1: Thermoelectric Refrigeration (T < T,) 


The optimum C.O.P. is given by the Carnot efficiency multiplied by 

a merit factor, Mc¢ which must be less than unity due to irreversible 
dissipation processes such as thermal conduction, 
r N (T/T) 


= = = X Me: Mc = - N+1 
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In addition (T> — T)max = $Zo7. 
Case 2: Thermoelectric Heating (T > To) 


Case 3: Thermoelectric Power Generation (T > or 


The maximum technical efficiency, no, occurs for a value of m = N 
and is given by the Carnot efficiency times a merit factor, M¢, which 
must be less than unity. Thus, 


To 
no = T Xx Me; Me= 


N-1 
N + T/T 


Among the most important features of these results are the following : 

The coefficient of performance (efficiency) of a thermoelectric heat 
pump increases with decreasing temperature difference between the 
junctions, while the efficiency of the generator increases with increasing 
temperature difference. In both the heat pump and generator, the 
efficiency depends on a merit factor, /, which expresses the fraction 
of the thermodynamic or Carnot efficiency actually achieved. The 
extent to which M is less than unity is a measure of the efficiency loss 
due to irreversible energy exchanges such as Joule heat (J?R) generation 
and thermal conduction. 

In both generators and heat pumps, if the cross section areas of the 
semi-conductor legs have been adjusted to give optimum performance, 
then their no or &) values depend on a single parameter expressing the 
properties of the materials. This factor is called the figure of merit Z, 
and is written Z, = Z = S’o/x, for a specific material, or a single leg of a 
thermocouple. In materials research for thermoelectric applications, it 
is thus Z which has to be maximized over the whole temperature range 
of application. 

The C.O.P. (& 9) and generator (no) efficiencies are independent of the 
volume of semiconductor in the thermoelement. Thus the thickness / 
of the semiconductor legs may be decreased without altering 6 or 70 
until the semiconductor-metal contact resistance becomes a significant 
source of energy loss from Joule heating. 

Some typical results of this phenomenological treatment for the 
temperature range covered by intermetallic semiconductor compounds 
are as shown in Tables I and II. 


3. SOME TECHNICAL ADVANTAGES OF THERMOELECTRIC HEAT PUMPS AND GENERATORS 


The prospects for widespread use of thermoelectric heat pumps are 
based essentially on certain characteristics that are summarized below : 
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TABLE I.—Thermoelectric Heat Pump Performance at Ty = 300° K. 


Coefficient of Performance of a single stage thermoelectric heat 
pump used as a refrigerator or as a heater. 7) = 300° K. 
AT (cooling) = 7) — 7, AT (heating) = T — T». 


Temperature C° 5 10 20 25 30 40 50 
difference F° 9 18 36 45 54 72 90 


Figure of merit, Zo = 5.0 & deg™! 


C.O.P. Cooling eA. 61 27 2.0 1.6 1.0 0.7 
Heating 13.6 7.3 4.0 a5 2.8 2.3 2.0 


Zo = 3.0 X 10-* deg 
C.O.P. (&) Cooling 8.9 4.1 1.8 1.3 
Heating 52. 3.0 25 22 1.8 1.6 
Zo = 2.5 X 10° deg™ 


Heating So 46 26 23 1.8 1.6 1s 


- 


Zo = 2.0 X 10°? deg 


C.O.P. Cooling 63 2.8 


Heating 7.0 4.0 3 2.0 4.7 1.5 1.25 


= 


Zo. = 1.5 deg 


C.O.P. (&) Cooling 49 2.2 083 0.55 0.33 0.12 tae 
Heating 1.8 1.4 1.1 


TABLE I1.—Thermoelectric Power Generation at Moderate Temperatures. 


Efficiency, no, of a single-stage generator. 7) = 300° K. 
T 
AT = 7 To, M at 
Temperature C° 100 200 300 400 500 700 
difference F° 180 360 540 720 900 1080 


Figure of merit, Zo = 5.0 10-3 deg 
a % 6.8 12.5 17.2 21.5 25.0 32.0 


Figure of merit, = 3.0 10°* deg”! 


no % 4.9 9.2 13.0 16.5 19.5 25.0 
Figure of merit, Z) = 2.0 deg”! 

«% 3.8 7A 10.0 13.0 15.5 20.0 
Figure of merit, Z) = 1.5 X 

no % 3.0 5.8 8.2 10.5 13.7 16.5 
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Thermoelements possess extreme technical simplicity. The ab- 
sence of motors and compressors results in a minimum of maintenance 
and the absence of lubrication problems, fluid working substance, and 
noise or vibration. They thus possess high technical reliability com- 
bined with ease of automatic control and continuous variation of heat 
pumping capacity by simple current regulation. 

Only a simple current reversal is required to reverse the direction of 
heat pumping to provide heating in cooling installations and vice versa. 
Thus convenient excursions can be made from room temperature in 
either direction with a single miniaturized thermoelement or a large 
battery of thermoelements. 

Thermoelements have high coefficients of performance for applica- 
tions requiring small temperature differences (Table I). 

Thermoelectric heat pumps have greater efficiency than compressor 
units for small installations. Their efficiency will become comparable 
with compressor units for larger installations when Z values above 
5 & 10-* are obtained. 

Considerable adaptability in design is possible with heat pumping 
thermoelements; thus, panels of thermobatteries or individual thermo- 
elements may be easily integrated with the most efficient heat exchange 
regions. 

The first order theory shows that the C.O.P. of the thermoelectric 
heat pump is independent of the volume of semiconductor material. 
Thus, only a small thickness of semiconductor may eventually be re- 
quired in thermoelements providing contact resistances can be suffi- 
ciently reduced. 

Thermoelectric heat pumps can work effectively under extremes of 
ambient temperature or in corrosive environments where mechanical 
units would not function easily or reliably, for example, arctic, tropical, 
outer space, or missile environments. 

Thermoelement fabrication should lend itself to inexpensive mass 
production techniques, since polycrystalline materials can be used with 
doping additions controlled to only moderate precision. 

The chemical elements involved in many suitable materials are often 
inexpensive and are readily obtainable in large quantities. 

An inherent disadvantage of thermoelectric heat pumps in household 
applications is the need for d.c. in their operation. An inexpensive 
germanium rectifying unit with low ripple content is thus necessary. 
On the other hand, the d.c. requirement is not a disadvantage in trans- 
port vehicles, or in laboratory and often factory installation, where d.c. 
is normally available. 

The advantageous features of thermoelectric generation are similar 
to those summarized for the thermoelectric heat pump. In the thermo- 
electric generator, however, the efficiency increases with increasing tem- 
perature difference. This characteristic often offsets the greater diff- 
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culty in obtaining high figures of merit with the presently investigated 
materials at temperatures in excess of 500° K. The reader is referred 
to Table II for the relative importance of these two factors. The most 
favorable characteristics of the thermoelectric generator follow, of 
course, from the elimination of the mechanical heat engine, and the 
direct and adaptable way in which thermoelements can be integrated 
into and around the heat source. Direct energy conversion is thus 
possible at the high temperature of the heat source, in contrast to the 
steam turbine and dynamo combination. Thermoelectric generation 
panels may be applied directly to large areas of a heated surface or 
thermoelements can be suitably integrated into pipes or baffles for heat 
exchange. 

The thermoelectric material used will of course depend on the tem- 
perature range to be utilized. Naturally, a thermoelectric generator 
may advantageously consist of two or even three different materials 
suitable in different temperature ranges. Different series connected 
materials, each maximized by proper doping additions, could form com- 
posite thermocouple legs. The different portions of each leg would 
then generate electricity at optimized efficiencies in different tempera- 
ture ranges, such as 2000 to 1200, 1200 to 700, and 700 to 300° K. 
As the temperature of the source, or combustion gas, decreased along 
the length of the heat exchange path, the higher temperature materials 
would be omitted in the thermocouples. Thus, such thermoelements 
would consist of successively three, two, and single-component semi- 
conductor legs. Materials used for the range below 700° K would be 
generally similar to those employed in thermoelectric heat pumps 
operating near 300° K. A figure of merit of 3.0 « 10-* deg gives a 
calculated generator efficiency (no) of 17 per cent with a temperature 
difference of 400° above 7) = 300° K. 

When discussing thermoelectric generation, the point is sometimes 
raised that each thermoelement develops only a very small voltage. 
However, it should be kept in mind, when comparing the thermoelectric 
generator and dynamo, that a single thermoelement, or a small group 
of thermoelements, is analogous in function to a single inductive turn 
on an electromagnetic generator. Thus a thermoelement, having a See- 
beck coefficient S=0.24 mV/° K., with T=800° K., and T,=300° K., 
generates 0.12 volt. If the thermoelement occupies an area of 2.0 cm?, 
then a thermoelectric panel of area 1.0 m* should generate 600 volts. 
With good heat exchange arrangements, a volume of only 1.0 m* con- 
taining thermoelement arrays may yield ten to twenty times this voltage 
in a series connected arrangement. When such an array of thermo- 
junctions is connected to produce maximum power in an external work- 
ing load, a generator efficiency (no) of 18 per cent is obtained for a figure 
of merit of 2.5 x 10-* deg-!, with a temperature difference of 500° K. 
from 300° K. 
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Admittedly the thermoelectric technical efficiency possible with the 
materials measured to date is not equal to that achievable by conven- 
tional means, except in quite small capacity applications. It should 
be stressed, however, that technical efficiency is, after all, only one of the 
factors which decide the more relevant end product efficiencies such as 
saleability, cost, weight-space, and convenience. Theremoelectric 
methods already present many novel and attractive features as well as 
new design possibilities which are not feasible with conventional tech- 
niques of heat pumping and power generation. This factor can fre- 
quently offset the former disadvantage. Naturally these characteristic 
advantages will make thermoelectric machines even more attractive for 
larger capacity applications if their efficiency equals or exceeds electro- 
mechanical engines. 

4. FIELDS OF APPLICATION 

The characteristic advantages of thermoelectric engines lead to the 

promise of applications in many fields of industrial technology. We 


shall summarize some of these fields in broad categories in each of which 
exists a large number of different applications. 


Heat Pumping 


(a) Domestic, catering, and industrial refrigeration. 

(b) Air conditioning and heating. 

(c) Water cooling, drinking fountains. 

(d) Numerous novel small household appliances involving small 
scale heating and cooling, for beverages, food and laundering. 

(e) Household and hot water heating. 

(f) Special cooling and heating applications to aeronautical, marine, 
missile and power reactor technology. 

(g) Numerous chemical, agricultural, medical, dental, and scientific 
instrument applications. 

(h) Thermostating and temperature control of electric components 
and circuits. 

(¢) Air conditioning, cooling, and heating applications in trans- 
portation. 

(7) New fields of heat pump application arising out of the unique 
features of thermoelements. 


Power Generation 


(a) Nuclear power plant generation. 

(b) Utilities and waste heat utilization. 

(c) Generation of electricity from natural gas. 

(d) Transportation : thermoelectrically-powered electric motor units 


May, 1959.] THE FRANKLIN INstITUTE LABORATORIES 451 


for propelling ships, submarines, locomotives, etc., using nuclear 
or hydrocarbon generated heat. 

Exhaust or waste heat power generation in transport vehicles 
(automobiles, tractors, locomotives, ships, submarines, aircraft, 
missiles), for battery charging, ignition, motor driving thermo- 
electric cooling, instrument and radio powering, etc. 

Solar generation. 

Portable generators of moderate power (less than 2000 watts) 
for isolated situations. 

Portable, very low power generators (less than 100 watts) for 
communication and other electronic applications. 


5. PRESENT SITUATION REGARDING FIGURES OF MERIT 


The results shown in Tables I and II stress the importance of 
fabricating materials of high figure of merit. In the low temperature 
range, —50 to 300° C., the best materials so far reported lie in the V—-VI 
and IV-—VI systems of semiconductor compounds. Solid solutions of 
the compounds within a system are often used to reduce the thermal 
conductivity of the materials. ‘Typical examples are suitably doped 
alloys of PbTe and PbSe(1) or Bi.Te; and Bi.Se;(2). Power generators 
having an efficiency between 8 and 10 per cent have also been made 
using these materials in the temperature range below 300°C. High 
temperature generation materials such as the mixed valence oxide com- 
pound semiconductors or systems of III-V compounds such as InAs- 
GaAs solid solutions show promise but have not yet been reported to 
have efficiencies greater than 10 per cent over temperature ranges ex- 
tending above 500° C. It is important to stress, however, that ma- 
terials research in this field is in its infancy. Only a small fraction of 
the large number of possibly useful materials has been investigated and 
in existing semiconductor compounds of promise such factors as oriented 
crystallization and defect-free crystals have been relatively neglected. 


6. CONCLUSION 


There seems little doubt that the impact of thermoelectric methods 
in refrigeration, heating, power generation, and transportation will be 
felt with increasing force as time proceeds. Obviously, the commercial 
realization of the developments in this field depends vitally on vigorous 
research and development in solid state materials. Moreover, research 
and development may be expected to be well rewarded in a field as 
young as this. It is thus not inappropriate to conclude this discussion 
by broadly summarizing the basic areas for unrelenting effort in research 
and engineering design, through which increasing improvement in per- 
formance can be attained. 
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Principal Areas for Research and Development 


The development of materials of increasing figure of merit over the 
whole temperature range of application. 

The development of techniques for further reducing semiconductor- 
metal contact resistances for miniaturized thermoelements. 

The arrangement of mass production methods for materials fabrica- 
tion and thermoelement formation. 

The design of heat exchange arrangements of increased technical, 
weight, and cost efficiency. 
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SOME PROBLEMS OF CHEMICAL KINETICS AND 
Reactivity. Vov. I, by N. N. Semenov, 
translated by J. E. S. Bradley. 305 pages, 
diagrams, 54 X 84 in. New York, Perga- 
mon Press, Inc., 1958. Price, $7.50. 


Nobel prize winner Semenov has combined 
an earlier, limited-edition, Russian book and 
a collection of U.S.S.R. Academy of Science 
papers in preparing this edition. It contains 
a very critical review of the latest interna- 
tional literature and much information from 
Russian studies that were previously unpub- 
lished. This is not intended as a complete 
treatise of the subject. It is replete with 
personal views, criticisms, and suggestions. 
Many areas requiring research are enumer- 
ated. The author urges broader understand- 
ing of Kinetics and reactions through basic 
research and general investigation rather than 
the study of isolated observations. ‘The book 
is divided into three parts. 

Part I deals with radical reactions, includ- 
ing chapters on uniradical reactivities, alter- 
nate univalent free-radical reactions, and reac- 
tions of biradicals. Topics included in the 
discussion are activation and bond energies, 
polar and steric factors, radical decomposition 
and isomerization, intermediate and chain re- 
actions, and an interesting treatment of meta- 
stable states in complex molecule photoreac- 
tions. The treatment of biradicals in the 
kinetics of chain reactions provides a thought- 
provoking excursion into this important sub- 
ject. 

Part II, Initiation and Termination of 
Chain Reactions, includes chapters on molec- 
ular dissociation and radical combination, ions 
of variable valency as chain initiators, and 
wall initiation and retardation of chain reac- 
tions. This reviewer found the last named 
chapter to be a highlight. The author ex- 
plores many facets of wall effects including 
free radical generation and removal, free va- 
lencies, radical production, concepts pertain- 
ing to heterogeneous catalysis, etc. Support- 


ing information of a quantitative nature has 
been combed from many sources and con- 
cisely presented. 


BOOK REVIEWS 


Part III, Kinetics of Chain Reactions, 
treats with chain reactions competing with 
reactions between saturated molecules. En- 
ergy dependence, relative rates, specific de- 
compositions and free-radical production are 
described. 

Appendices treat with the activated com- 
plex and the quantum-mechanical calculation 
of activation energies. There are many refer- 
ences for every chapter and these are supple- 
mented by a name and subject index. 

Space limitations preclude a more detailed 
review. Suffice to say that scientific person- 
nel concerned with kinetics or reaction mech- 
anisms will find this a highly worthwhile, 
stimulating contribution. 

RoBertT A. BAKER 
The Franklin Institute Laboratories 


HUMAN ENGINEERING, by Ernest J. McCor- 
mick. 467 pages, illustrations, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 
1958. Price, $8.00. 


This text has been directed toward the con- 
sumers of human engineering recommenda- 
tions—design engineers and project managers. 
Rather than present an oversimplified catalog 
of those human characteristics which are 
important when men act in concert with ma- 
chines, the author has chosen to give the 
reader a certain amount of insight into the 
process by which these human characteristics 
were determined. This background material 
has been presented clearly and adequately. 
Psychological nomenclature of parochial ap- 
peal has been kept down, and recommenda- 
tions are presented in a well reasoned manner. 

As with any text of this nature, there are 
omissions and weaknesses, although the au- 
thor has done a fine job of culling the unclassi- 
fied literature of government reports. Never- 
theless, there are strong trends in current 
human engineering thinking which have not 
been treated. One such trend is the ever 
increasing recognition of the complementarity 
of human engineering and operations research 
in the study of large and complex manned 
This concept deserves a detailed 
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and extended coverage. By and large, how- 

ever, McCormick has prepared a useful, well 
written book on an infant technology. 

Ezra KRENDEL 

The Franklin Institute Laboratories 


ENGINEERING VIBRATIONS, by Lydik Jacobsen 
and Robert S. Ayre. 564 pages, diagrams, 
6 X9in. New York, McGraw-Hill Book 
Co., Inc., 1958. Price, $10.00. 


This is a welcome and useful treatise on 
vibration analysis, particularly adaptable for 
directing the solution of a large class of typical 
engineering vibration problems. 

The work is designed as a suitable course 
on an early graduate level and includes an 
extensive bibliography. The unique feature 
of the book is a new approach in the setting 
up of the differential equations and actually 
carrying out the solution for an important 
class of engineering cases. By this method, 
the authors have successfully developed a true 
physical feel of the subject. Along with the 
analysis, free use of graphical and tabulated 
methods are used, such as time displacement 
curves for damping ratios, and the establish- 
ment of roots in transcendental frequency 
equations, as well as the picturing of modes 
of vibrations of multi-degree-of-freedom sys- 
tems. 

The early chapters cover linear and non- 
linear single-degree-of-freedom systems, tran- 
sient and steady forced vibrations with and 
without damping, including approximate solu- 
tions. 

Part II covers linear undamped multi- 
degree-of-freedom systems with various types 
of couplings. Numerical methods, with ac- 
tual solutions for multi-mass systems, include 
the Holzer and the Stodola-Rayleigh methods. 
The shearing vibration of a six story building 
is discussed. 

The last chapters deal with continuous sys- 
tems, as the stretched string, bars with various 
dampings, pressure waves, and such matters 
as the vibration of shear beams. 

The authors are well known in their con- 
tributions to vibration analysis, and in par- 
ticular to vibration problems in structures, 
and this work reflects a senior approach to 
the subject. They are to be complimented 
in presenting a difficult subject in a clear-cut 
simple manner. 

RUPEN EKSERGIAN 
The Franklin Institute Laboratories 
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PRINCIPLES AND APPLICATIONS OF RANDOM 
NotsE THEORY, by Julius S. Bendat. 431 
New York, 


pages, diagrams, 6 X 9 in. 
Price, 


John Wiley & Sons, Inc., 1958. 

$11.00. 

The sciences of communication and control 
have gone through a complete evolution cycle 
from a primitive status, characterized by a 
purely deterministic approach, to the recent 
developments resulting from the recognition 
of the statistical nature of the processes in- 
volved. 

One of the results of this evolution is the 
appearance of books intended to recapitulate 
and clarify the major classical original con- 
tributions illustrating at the same time the 
various applications. 

The difficulty of locating and interpreting 
the numerous and not always coordinated 
original works makes the availability of well 
written books of this kind a blessing to scien- 
tists, engineers and graduate students. 

Dr. Bendat’s work is especially devoted to 
the study of the principles and applications 
of random processes rather than to that spe- 
cific body of knowledge which goes under 
the name of “Information Theory". The re- 
viewer agrees with the author’s opinion that 
this distinction is expedient although regret- 
table in view of the areas of overlap between 
the two philosophies. 

The book is divided into ten chapters and 
it is pleasing to find detailed summaries of 
the individual chapters at the beginning of 
the book in addition to the conventional index. 
Furthermore six chapters end with a conclud- 
ing summary that is highly effective and clar- 
ifying. Appendices are also found at appro- 
priately chosen points. 

The first chapter gives an introduction to 
random processes, correlation functions and 
system analysis. The concepts of stationarity 
and ergodicity are presented in an elementary 
way yet with sufficient sophistication to warn 
the reader against indiscriminate acceptance 
of certain hypotheses. Systems are classified 
in terms of linearity, time-invariance of pa- 
rameters, finiteness of operating time. A 
very brief treatment of superposition in the 
time domain and of the related frequency 
response functions is also given. 

Chapter 2 deals with power spectra and 
correlation functions as time averages. The 
presentation of this material proceeds from 
the periodic to the aperiodic and, eventually, 
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to the random case. The sampling theorems 
are presented with special attention to their 
relation to the resolving power in time and 
frequency for any physical system. The 
use of correlation techniques to extract ‘‘hid- 
den” periodicities in quasi-random processes 
and to obtain input-output relations for linear 
constant parameter systems is discussed. 

Chapter 3 opens with an analysis of en- 
semble averages followed by a condensed ex- 
position of probability theory culminating in 
the Central Limit Theorem for one random 
variable The theory is extended to multi- 
variable case with the presentation of distribu- 
tion of sums, characteristic and moment gen- 
erating functions. The effects of a change 
of variable in probability density functions 
is thus discussed and the chapter ends with 
the presentation of the multivariate normal 
distribution, the Central Limit Theorem in 
dimensions and numerous applications to the 
analysis of random noise. 

Chapter 4 is devoted to Optimum Linear 
Prediction and Filtering. Although the least- 
mean-square error philosophy is given due 
predominance a special effort is made to com- 
pare the various conditions and assumptions 
chosen by different authors on account of 
mathematical expediency or physical signifi- 
cance. (See Table 4-1 and Section 3 of Chap- 
ter 4). The historical approach together with 
frequent examples will help the reader to reach 
a firm understanding of the theory for any 
combinations of: constant, or time-varying, 
parameter linear systems, finite, or infinite, 
operating-time systems and stationary, or 
non-stationary, random processes. 
attention is also directed to 
and the differences between 


The reader’s 
the similarities 
the filter-predictor problem and that of the 
extraction of a signal from noise under as- 
sumptions that may differ from the least- 
mean-square error criterion. 

Chapter 5 analyzes theexponentially damped 
autocorrelation function which has proven 
useful in many cases. Examples from ther- 
mal, tube and circuit noise, telegraph and 
other signals, fading phenomena in turbulent 
media are cited to support the above point 
of view. It is shown that a suitable combina- 
tion of several exponential, or exponential- 
cosine, functions can be used to fit many 
physical cases. A discussion of Markoff proc- 
esses closes the chapter, 
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Chapter 6 is devoted to analog computers 
designed for driving forces characterized by 
exponentially damped cosine functions and 
for random errors in the initial conditions 
that are assumed to be statistically independ- 
ent among themselves and from the driving 
forces. 

Chapter 7 offers quantitative criteria to 
estimate the statistical errors occurring in the 
measurement of correlation functions with 
finite bandwidth and finite integration time. 
The effect of a varying correlation interval 
and the concepts of “no correlation’ and 
“false correlation”’ probabilities are presented. 

Chapter 8 is devoted entirely to the problem 
of the demodulation from RF to video of 
signals in noise. Correlation detectors are 
compared to envelope (‘‘linear’’) and square- 
law detectors. The effects of automatic gain 
control are briefly mentioned. The impor- 
tant role of a finite bandwidth-integration 
time product is again stressed in this Chapter. 

Chapter 9 formulates the time-variable pa- 
rameter filter problem in the sense of opti- 
mizing, at the filter output, an appropriate 
signal function or some transformation thereof 
when the system input is a linear combination 
of signals and noise. The adopted optimiza- 
tion criterion is to minimize the mean-square 
ensemble difference between the ideal and the 
effective outputs. This difficult problem is 
handled both from an analysis and a synthesis 
point of view and a general solution is ob- 
tained. 

Chapter 10 is devoted to the so called zero- 
crossing problem which is important in many 
applications. The usefulness of the expected 
number of zero crossings and of the output 
power spectral density is considered in detail 
with reference to pure random-noise inputs 
or to combinations of signals and noise. 

A bibliography of 141 references and an 
alphabetical index close the book 

The well coordinated presentation of known 
material, the original contributions, and, last 
but not least, the author's successful effort 
to organize the book in a way which should 
facilitate the reader's task of absorbing a large 
amount of knowledge, deserves the highest 
praise especially in an epoch like ours in which 
cryptic writing is often confused with virtu- 
osity. 

P. L. BARGELLINI 
University of Pennsylvania 
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A History OF SCIENCE, by George Sarton. 
554 pages, diagrams, 64 X 9} in. Cam- 
bridge, Harvard University Press, 1959. 
Price, $11.00. 


This work, completed just before Dr. Sar- 
ton's death, is his last book. Although the 
subtitle would indicate that the volume deals 
only with Greek science in the last three cen- 
turies B.C., it does in fact cover Roman culture 
and Latin letters, as well as the culture of 
Eastern Europe, Egypt, and Western Asia. 
The three centuries treated witnessed the 
establishment of geometry, astronomy, anat- 
omy, and grammar, and the beginnings of 
technology and medicine. Dr. Sarton evalu- 
ates each ancient achievement in relation to 
what bearing it had on later achievements (in 
contrast to his original plan). 

The book is divided into two parts. Part 
One treats the third century B.C., including 
the Alexandrian renaissance, Euclid, the be- 


ginnings of astronomy, Archimedes and the 
mechanics, geography, physics and technology 
(public works), anatomy, medicine, the li- 
brary, philosophy and religion, history and 


language. Part Two, on the last two cen- 
turies B.c., deals with the growth of Rome and 
continues the study of the subjects enumer- 
ated for Part One. 

This monumental work will long stand as 
the authoritative reference work in the field. 
It is a matter of deep regret that Dr. Sarton’s 
original plan—to cover the history of science 
through the 19th century—cannot now be 
carried out. 


RANDOM VIBRATION, edited by Stephen H. 
Crandall. v.p.,diagrams,6 X 9in. Cam- 
bridge, The Technology Press, 1958. Price, 
$10.00 (paper). 

This volume was prepared for use in the 
special summer program on Random Vibra- 
tion, offered at M.I.T. in 1958. The problem 
of random vibration is of vital importance 
in the high power rocket and jet engine fields, 
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since it severely affects not only the perform- 
ance of the vehicle, but also that of the air- 
borne electronic equipment. The several 
papers included in the volume describe the 
new concepts required to extend ordinary vi- 
bration theory into the field of random vibra- 
tion. The authors also set forth the current 
state of the art of designing and testing equip- 
ment to withstand random vibration. The 
book is divided roughly into two parts: the 
first six chapters deal with basic concepts and 
background material; the last six chapters 
cover measurement of random vibration, de- 
sign problems in aircraft and missiles subject 
to random vibration, excitation from jet en- 
gines, rocket motor excitation, design of simu- 
lation equipment, and mechanical loading. 
This book, the first of its kind, should be 
required reading for all designers involved in 
rocketry, missiles and jet engine aircraft. 


Cosmic ELECTRODYNAMICS, by J. W. Dungey. 
183 pages, diagrams, 54 X 84 in. New 
York, Cambridge University Press, Inc., 
1958. Price, $6.00. 


One of the Cambridge Monographs on Me- 
chanics and Applied Mathematics, this small 
volume treats the relatively new subject of 
cosmic electrodynamics from a_ theoretical 
viewpoint. The author discusses the funda- 
mental physical theories involved and explains 
the mathematical methods used. After a pre- 
liminary discussion, the author covers these 
eight general subjects: the velocity-distribu- 
tion method; the motion of magnetic fields; 
magnetostatics; magnetodynamics; accelera- 
tion of charged particles to high energy ; solar 
phenomena (sunspots, corona, prominences, 
flares M regions) ;magneticstormsandaurorae ; 
and ionospheric electrodynamics. The mono- 
graph will be appreciated by those interested 
in cosmogony, astronomy, astrophysics, solar 
physics, atmospheric physics and others in 
the fields of cosmic ray work, the Earth’s 
magnetism and particle accelerators. 
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REFLEX Kiystrons, by J. J. Hamilton. 260 
pages, illustrations, 5} X 8}in. New York, 
The Macmillan Co., 1959. Price, $9.00. 


Written on an undergraduate level, this sur- 
vey incorporates an elementary outline of the 
theory and operation of reflex klystrons, in- 
cluding the fundamentals of cavity resonators, 
together with a discussion of the techniques 
of the design, manufacture and testing of such 
tubes. ‘Typical reflex klystrons, from some of 
the earliest production types to modern milli- 
meter wave tubes are described and illustrated. 
The last two chapters cover some unconven- 
tional models and the direction in which pres- 
ent developments are trending. 


FINE PARTICLE MEASUREMENT, by Clyde Orr, 
Jr. and J. M. DallaValle. 353 pages, dia- 
grams, 53 X 8}.in. New York, The Mac- 
millan Co., 1959. Price, $10.50. 

The authors have gathered in one volume 
the data on important techniques used in 
research laboratories for fine particle measure- 
ment. Covering theory and application, the 
book describes over seventy techniques cur- 
rently in use for determination of size, surface 
and pore volume. Ten chapters deal with 


various phases of fine particle measurement, 
including microscopy and sieving, sedimenta- 


tion, inertial techniques, radiation scattering 
and transmission, permeametry, gas adsorp- 
tion, liquid-phase sorption, and X-ray evalua- 
tion of pore size. 

Since fine particle measurement plays so 
important a role in widely diverse fields (for 
example, concrete, paints, foods, catalysis), 
the book should be of value to a wide audience. 
The presentation is limited to basic principles, 
use of apparatus, methods of data analysis, 
applications and limitations; it does not at- 
tempt to relate measurement with properties. 


THE TECHNICAL WRITER, an aid to the presen- 
tation and production of technical litera- 
ture, edited by J. W. Godfrey and G. Parr. 
340 pages, illustrations, plates, 5} X 8} in. 
New York, John Wiley & Sons, Inc., 1959. 
Price, $8.50. 


New information in all branches of science 
and technology needs to be disseminated as 
rapidly as developments take place. The de- 
signers and researchers often have not the 
time or perhaps the ability to write up prog- 
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ress in a way that could help others in indus- 
try. This creates a need for skilled technical 
writers who have that ability. 

The term “technical literature’ covers a 
variety of publications which have the com- 
mon object of conveying information as ef- 
ficiently and clearly as possible, and this is 
the underlying art of technical writing. In 
larger organizations the technical writer often 
forms part of a team producing technical liter- 
ature In smaller firms the engineer or de- 
signer may do his own writing and production 
To give the technical writer an understanding 
of the technique of presentation and produc- 
tion is the purpose of this book. The au- 
thors have assumed that the writer is already 
grounded in the principles of composition, and 
mainly needs to avoid common pitfalls and 
acquire a readable style. The book also 
stresses that the writer's responsibility for his 
work carries through the entire publication 
process. The book familiarizes him with the 
processes through which his manuscript passes, 
aiming to promote efficient, neat and econom- 
ical production. 

The examples are to a large extent drawn 
from electrical and electronic engineering prac- 
tice but the principles hold for all branches of 
science. The attention given here to circuit 
diagrams stems from the fact that in modern 
engineering, mechanical and electrical tech- 
nologies are interdependent, and nearly all 
allied sciences make use of electronic circuits. 

In its analysis of the forms of technical 
literature, the book describes the technical re- 
port, paper, article, leaflet, sales information 
and books. Editorial procedure, typography 
for technical publications, set-up of a technical 
publications unit, and technical authorship 
are each analyzed in detail. 

Appendices are added covering casting-off 
tables, proof-correction marks, printing terms, 
and associations for technical writers. 


THE ASTRONOMER’'S UNIVERSE, by Bart J. 
Bok. 107 pages, illustrations, 54 84 in. 
New York, Cambridge University Press, 
1959. Price, $3.75. 


Written in non-technical terms, this text 
gives an account of recent developments in 
physics, mathematics and radio astronomy 
that have put new tools and methods into the 
hands of astronomers. Selected topics include 
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the measurement of size, distance and bright- 
ness of stars, problems of the solar system 
and the Milky Way, radioactivity and map- 
ping the heavens. In the final chapter the 
author enters into the realm of theory and 


speculations. 


ANIMAL Lire, by Lorus J. Milne and Margery 
Milne. 367 pages, diagrams, illustrations, 
7 X 93in. New York, Prentice-Hall, Inc., 
1959. Price, $6.95. 


New vitality is given to the age-old science 
of life’s origin in this richly pictorial text. 
It is written with great clarity in a refreshing 
style beginning with the basic information 
The impor- 
and 


concerning the living processes. 
tance of coordination, both chemical 
nervous, in active multicellular animals is fol- 
lowed by contrasts of single-celled life. The 
dramatis personae of zoology is outlined for 


their economic importance. 


CONSTITUTIONAL DIAGRAMS OF URANIUM AND 
THor1uM ALLoys, by Frank A. Rough and 
Arthur A. Bauer. 153 pages, diagrams, 
6 X 9in. Reading Mass., Addison-Wesley 
Publishing Co. Inc., 1959. Price, $5.00. 


This valuable compilation of United States 
and United Kingdom uranium and thorium 
constitutional diagrams is a compact and 
highly useful source of reference material for 
workers in the field. ‘The book is divided 
into two major sections, uranium and thorium 
alloys, each preceded by a discussion of the 
transformation and melting temperatures of 
thebase metal. The varioussystemsare listed 
in alphabetical order, including both binary 


and ternary systems. 


ELEMENTS OF PuysicAL METALLURGY, by Al- 
bert G. Guy. 528 pages, diagrams, illustra- 
tions, 6 X 9in. Reading, Mass., Addison- 
Wesley Publishing Co., Inc., 1959. Price, 
$9.50. 


This second edition of a widely used text- 
book for introductory courses in physical met- 


allurgy contains several revisions in the intro- 


duction of dislocation theory. There is addi- 
tional material on electron theory of metals, 
in view of the increased importance of electric 
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phenomena in semiconductors. A new chap- 
ter is included on ‘‘Phases in Metal Systems” 
and other material designed to bring the text 
completely up to date. Written especially 
for science and engineering or metallurgy stu- 
dents since it treats basic principles. 


FILLER METALS FOR JOINING, by Orville T. 
244 pages, diagrams, illustra- 
New York, Reinhold 
Price, $7.00. 


Barnett. 
tions, 6 X 9 in. 
Publishing Corp., 1959. 


This volume summarizes the industrial prac- 
tice in selecting the proper filler metal for 
any metal joining situation. It covers all of 
the ferrous and non-ferrous joining metals. 
Welding, brazing and soldering are treated 
with major emphasis on products for manual 
and automatic arc welding. Several chapters 
are devoted to arc welding electrodes and 
automatic welding of mild and low alloy steels 
and discusses materials in an area where no 
specifications exist as yet. The chapters on 
new iron powder electrodes, the latest develop- 
ments in non-ferrous metals, surfacing and 
tungsten electrodes contain much material 
that has never been published before. An 
excellent reference for welding engineers and 
designers, 


OCEANOGRAPHY AND MARINE Brooecy, by H. 
218 pages, diagrams, illustrations, 

New York, The Macmillan 
Price, $7.00. 


Barnes. 
54 X 8} in. 
Co., 1959. 


In this volume the author comprehensively 
describes the ingenious instruments and meth- 
ods by which we obtain our knowledge of the 
sea. Methods of taking bottom cores and of 
using nets, pumps, trawls, dredges and grabs 
are all carefully explained. A section is de- 
voted to the properties of water itself— 
temperature, salinity and currents. Particu- 
lar attention is given to new techniques as for 
instance, the bathythermograph. The last 
chapters deal with the most spectacular of 
modern developments—underwater photogra- 
phy and television. Aerial photography is 
also included as a growing contribution to 
oceanography. Students of marine biology 
will find this volume of the greatest interest 
while marine research workers will value the 
assemblance of their techniques. 
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WELDING oF Ptastics, J. Alex Neumann. 
279 pages, diagrams, illustrations, 6 X 9 in. 
New York, Reinhold Publishing Corp., 
1959. Price, $7.25. 

Various physical and chemical properties of 
individual plastics presently used for welded 
construction are completely covered in this 
volume, with regard to their use in chemical 
and allied processing. All presently known 
techniques of welding are covered in detail 
and there are many step-by-step procedures 
with photographs. One of the most complete 
corrosion resistance tables yet published is 
included, plus a concise guide for the choice 
of an appropriate plastic construction material. 
This text will be of substantial help to engi- 
neers and manufacturers who are interested 
in present and potential applications in the 
process industries. 


LUBRICATION SCIENCE AND TECHNOLOGY, ed- 
ited by John Boyd. 346 pages, diagrams, 
illustrations, 8} X 11} in. New York, Per- 
gamon Press Inc., 1959. Price, $10.00. 


This book represents Volume II in the 
Transactions of the American Society of Lu- 
brication Engineers, constituting a valuable 
group of papers on widely diverse fields of 


lubrication. The papers are all written by 
leaders in their respective fields giving a source 
of reference not only for those concerned di- 
rectly with lubrication, but also to many in- 
dustrial research engineers and designers. 


THE ATOM AND THE ENERGY REVOLUTION, 
by Norman Lansdell. 200 pages, plates, 
6 X 94 in. New York, Philosophical Li- 
brary, 1958. Price, $6.00. 

In this book Lansdell presents a balanced 
study of the role of atomic energy in the world 
of today and tomorrow. He concerns him- 
self only with the peaceful uses of energy, 
including the atom. ‘The issue is the reorien- 
tation of human life and society, the impact 
of power development on other spheres of 
activity (on other industrial techniques, pri- 
vate industry, and the balance of world trade). 
The book is addressed to the specialist to ac- 
quaint him with other disciplines, to the busi- 
ness man who must adjust to the new indus- 
trial and commercial factors introduced by 
atomic energy, to the citizen to show him 
the social and political implications of the 
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new sources of energy, their potentialities, 
demands, and possible cost. 

The serious problems which must arise dur- 
ing the next fifty years, arising from the world 
pattern of energy demand, are neither the 
same nor equally crucial all over the world. 
New sources of energy are required on dif- 
ferent scales (large stations delivering power 
and/or heat, medium and small fixed units, 
and mobile units). Such sources are wind and 
water power, vegetable resources, solar en- 
ergy, electro-chemical generation of energy, 
direct conversion of radiation into electric 
current, nuclear fission, nuclear fusion, and 
the annihilation of matter. Present atomic 
energy developments relate mainly to large 
units. 

The structure of the atom and its energy 
are analyzed. The study of the methods of 
releasing atomic energy is illustrated with 32 
plates. 

Lansdell shows the sources of natural ma- 
terials for atomic energy development, the 
exploitation of atomic energy, the political and 
commercial organizations for atomic energy 
development, national and international, radi- 
ation risks, and insurance against them. 


OF OPERATIONS 
Saaty. 421 
New York, 
1959. Price, 


MATHEMATICAL METHODS 
RESEARCH, by Thomas L. 
pages, diagrams, 6 X 9 in. 
McGraw-Hill Book Co., Inc. 
$10.00. 


the scientific approach 
greater 


Operations research 
to operational problems for the 
fulfillment of objectives—is a parallel develop- 
ment to an increased realization of the use- 
fulness of staff planning and analysis func- 
tions. The implementation of these functions 
has led to rapid development of new meth- 
odology, and application of old methodology 
tonew problems. The book requires for back- 
ground a course in calculus, some advanced 
calculus and rudimentary matrix theory, and 
some analytical experience. 

Part I of this study gives background per- 
spective. The chapter on mathematical mod- 
els illustrates, by brief discussions and ex- 
amples, how the need for such tools arises, and 
some types of related theoretical questions. 

Part II treats of optimism and some of its 
problems, linear and quadratic programming, 
and the theory of games. 
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Part III covers probability, applications, 
fundamental statistics and their problems, and 
a resume of the queueing theory. There are 
many references for further information on 
the subject. 

Part IV is an essay covering some thoughts 
on creativity and elementary problems and 
projects. In the final analysis, the individual 
himself, his talents and methods of approach- 
ing problems, form the fundamental qualities 
required in solving problems. The analytical 
methods studied give a reliable approach to 
problems. Many of man’s problems have 
been successfully solved by a combination of 
intuitive talent and a desire to accept the 
challenge presented by a new problem. There 
is a place for both the analytical and the 
intuitive methods of attack. The mathema- 
tical methods are given primarily to increase 
understanding of problems and secondarily to 
give solutions to these problems. 


Basic GEOLOGY FOR SCIENCE AND ENGINEER- 
ING, by Edward C. Dapples. 609 pages, 
illustrations, 6 X 9 in. New York, John 
Wiley & Sons, Inc., 1959. Price, $9.50. 


This book presents a systematic analysis of 
the processes which are considered basic to 
the understanding of physical geology. It ex- 
plains geologic processes, and the techniques 
of geologic interpretation and reconstruction 
of former geologic conditions. 

The method of study of physical geology 
here is (1) observation and grouping of data, 
(2) classification of data into forms which can 
be expressed in graphs and tables to determine 
any measurable relationship between proper- 
ties, (3) interpretation of the relationships to 
form significant generalizations of the behavior 
of geologic processes. The emphasis is on the 
application of geology rather than the descrip- 
tion of localities where specific conditions are 
known to occur. Dapples believes that there 
are certain aspects of physical geology which 
must be considered basic to further study in 
the applied, descriptive or theoretical form. 
The pattern of his presentation is organized 
around empirical relationships. 

The body of scientific knowledge ranges 
from the taxonomy of fossil organisms to theo- 
retic analysis. The difficulty in the presenta- 
tion of the subject matter of geology to the 
scientist and the engineer is the selection of 
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what is to be considered the essential ingre- 
dients of the science. 

Dapples studies soil materials, physical and 
chemical properties of rock materials, rock- 
forming minerals, igneous rocks, soil-forming 
processes, streams, shoreline processes, ground 
water, wind deposits, processes associated with 
glaciation, sedimentary rocks, crustal defor- 
mation, and metamorphism. The study is 
enhanced with tables for the determination 
of common minerals. 


PLastic ANALYSIS OF STRUCTURES, by Philip 
G. Hodge, Jr. 364 pages, diagrams, 6 X 9 
in. New York, McGraw-Hill Book Co., 
Inc., 1959. Price, $10.50. 


For many generations, engineers have based 
the analysis and design of structures on a 
lineal theory of elasticity. Hodge concerns 
himself with plastic analysis to point out some 
of the disadvantages of elastic-design methods. 
Ductile structural materials such as steel and 
aluminum can withstand strains much larger 
than those encountered within the elastic lim- 
it, and carry stresses above the yield stress 
(strain hardening). This ductility in redun- 
dant structures permits a redistribution of 
stresses beyond the elastic limit, enabling the 
carrying of very considerable additional loads. 
Although the average stress may be below 
the elastic limit, highly stressed local regions 
may become plastic, making elastic analysis 
frequently unrealistic. 

Plastic design partially answers the objec- 
tions to elasticity. It takes full advantage 
of ductility. Nonuniform stress distributions 
can be advantageously redistributed. The 
method is mathematically simple. Through 
his analysis in this book, Hodge shows that 
these claims are fully justified in many prac- 
tical situations, and he recognizes that there 
are situations where plastic analysis does not 
fit. He explains how to apply plastic methods 
but also judges when these methods should be 
applied. 

To simplify presentation, Hodge has di- 
vided the book in two parts. Part | on beams 
and frames treats of the application of plastic 
methods to frame-type structures whose prin- 
cipal strength lies in their resistance to bend- 
ing. In addition to a strictly plastic or limit 
analysis, he deals with topics such as elastic 
plastic deformations, variable and repeated 
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loading, and direct procedures of design. The 
discussion of framed structures, where tensile 
stresses produce both an axial force and a 
bending moment, makes the transition to the 
combined stress problems of part II. This 
includes combined stresses in beams, circular 
plates, and circular cylindrical shells. Less 
complete treatment is given to more general 
plate and shell problems and to problems in 
plane stress. In conclusion, there is an intro- 
duction to some of the problems encountered 
in the dynamic loading of plastic structures. 

Asa teacher, Hodge has written this book to 
appeal both to the undergraduate or the first- 
year graduate student learning about struc- 
tures for the first time, and to the practicing 
engineer who has completed his course work 
years before. With this double aim he sepa- 
rated the theoretical material, such as proofs 
and theorems, into appendices. The second 
part, concerned with more complex structures, 
is written at a more advanced level. He sug- 
gests additional references where topics have 
not been fully covered. Problems are offered 
at the end of each chapter to encourage the 
reader to seek out original sources. 


LeGa Aspects OF CONSTRUCTION, by Walter 
C. Sadler. 387 pages, diagrams, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 
1959. Price, $8.50. 


This book, written for practical use in the 
construction industry, is focused primarily on 
a discussion of the rights and liabilities on 
construction, as a guide to help avoid legal en- 
tanglements and litigation in the areas of pos- 
sible conflict of interest. The development 
of the book has been directed to the business 
of contracting, with its many entwining rela- 
tionships, the relative liabilities of the many 
parties in the construction industry, and the 
property rights of abutting or adjacent parcels 
of land incidental to construction. The pri- 
mary parties are four-fold : owners, financiers, 
designers, and contractors. 
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The business of contracting is traced through 
licensing of contractors, mechanic's liens, in- 
surance, labor and management, the union 
contract, workmen’s compensation, compen- 
sable accidents, and arbitration. 

Further help in avoiding conflict of interest 
is given by showing court decisions in typical 
cases in such areas as architectural practice, 
architect-owner relationships, owner’s respon- 
sibilities, defective design and construction, 
general liabilities, and professional registra- 
tion. 

To avoid legal entanglements ‘‘at the bound- 
ary,’’ Sadler deals fully with boundary rights, 
structures at the boundary, lateral support, 
foundations, statutes and ordinances, money 
damages, and injunction and mandamus. 

Sadler, being an attorney and an engineer- 
ing and construction consultant, offers this 
book as a guide to ‘‘preventive law” for the 
construction industry. 


PuysicaL Laws AND EFFEcts, by C. F. Hix, 
Jr. and R. P. Alley. 291 pages, diagrams, 
5¢ X 9in. New York, John Wiley & Sons, 
Inc., 1958. Price, $7.95. 


This unusual book has been prepared to 
provide a centralized source of information 
on physical laws and effects. The authors, 
realizing that new ideas are frequently in- 
spired by old laws, have gathered new laws 
and old ones, often getting their material from 
half-forgotten sources. In a brief introduc- 
tion, the authors explain why such a compila- 
tion is useful. An 8-page index of physical 
Jaws and effects follows this introduction, with 
page numbers keyed to each, referring to the 
descriptions given in Chapter 3. Chapters 
4 and 5 are cross references, respectively, to 
physical quantities and fields of science, for 
the material in Chapter 3. A final chapter 
contains several rules, effects and laws which, 
because the authors do not vouch for their 
accuracy, are not included in the main body 
of the text. 


ag 
: 
: 
i 


462 Book Notes 


PUBLICATIONS RECEIVED 


ELECTRONIC FUNDAMENTALS AND APPLICATIONS, by John D. Ryder. Second edition, 721 
pages, diagrams, 6 X 9 in. New York, Prentice-Hall, Inc., 1959. Price, $13.35. 


PLANE TRIGONOMETRY, by Gordon Fuller. Second edition, 281 pages, diagrams, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 1959. Price, $4.75. 


FUNDAMENTALS OF ENGINEERING DRAWING, by Warren J. Luzadder. Fourth edition, thor- 
oughly revised, 720 pages, diagrams, 6 X 9 in. New York, Prentice-Hall, Inc., 1959. 
Price, $10.00. 


GENERAL Cuemistry, by William H. Nebergall and Frederic C. Schmidt. Alternate edition 
of “General Chemistry”’ by the same authors, 723 pages, illustrations, 7 X 9 in. Boston, 
D. C. Heath and Company, 1959. Price, $7.25. 


ELECTRICAL ENGINEERING: THEORY AND Practice, by William H. Erickson and Nelson H. 
Bryant. Second edition, 614 pages, diagrams, 6 X 9 in. New York, John Wiley & Sons, 
Inc., 1959 Price, $8.00. 


Mecuanics. Part II: Dynamics, by J. L. Meriam. Second edition, 420 pages, diagrams, 
6X9 in. New York, John Wiley & Sons, Inc., 1959. Price, $5.00. 


SEMICONDUCTOR ABSTRACTS, compiled by Battelle Memorial Institute. Volume I1V—1956 
issue, edited by E. Paskell. 456 pages, 8} X 11 in. New York, John Wiley & Sons, Inc., 
1959. Price, $12.00. 


Hi-Fi Furniture, by Jeff Markell. 224 pages, illustrations, 54 X 84 in. New York, Gerns- 
back Library, Inc., 1959. Price, $2.90 (paper) ; $5.00 (hard cover). 


INTRODUCTION TO SYMBOLIC LoGic AND Its APPLICATIONS, by Rudolf Carnap. Translation 
of “Einfiihrung in die symbolische Logik,’”’ revised and published for the first time in 
English. 241 pages, 54 X 8 in. New York, Dover Publications, Inc., 1958. Price, 
$1.85 (paper). 


Time's Arrow AND Evo.ution, by Harold F. Blum. Second edition, 224 pages, 6 X 9 in. 
Princeton, Princeton University Press, 1959. Price, $1.75 (paper). 

ORGANIC SYNTHESES. VOLUME 38. 1958. 120 pages, diagrams, 53 X 9 in. New York, 
John Wiley & Sons, Inc., 1958. Price, $4.00. 
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Durable Pleating and Creasing Proc- 
ess for Wool.—Wool trousers and 
skirts with creases and pleats that 
never lose their press will soon be avail- 
able all over America. Development 
work on volume production techniques 
for applying the new Si-Ro-Set process 
for the durable pleating and creasing 
of woolens and worsteds has now 
been completed, according to Max 
F. Schmitt, president of The Wool 
Bureau, Inc. 

Experimental work on the adapta- 
tion of the process to U. S. factory 
production has been going on in the 
United States since early 1958, shortly 
after its development in Australia. 
Si-Ro-Set was the result of research 
conducted by the Commonwealth Sci- 
entific and Industrial Research Or- 
ganization in that country. 

Mr. Schmitt said that two types of 
test production have been carried on 
under the supervision of Dr. Gerald 
Laxer, The Wool Bureau’s research 
director. 

One method involves the manual 
operation of conventional spraying and 
pressing equipment specially adapted 
for use with the new process. The 
other method utilizes new, automatic 
spraying and pressing equipment de- 
veloped specially for the Si-Ro-Set 
process under Dr. Laxer’s direction. 

Full information on how to use Si- 
Ro-Set was disclosed in March, in a 
technical booklet prepared for release 
to the industry. The Wool Bureau 
has also produced a motion picture to 
serve as a graphic method of instruc- 
tion in the application of the process. 

Personal instruction for clothing 
manufacturers is available through 


and Charles Macaluso 
at The Wool Bureau. Mr. Macaluso 
was added to the technical staff 
of The Wool Bureau at the end of 
1958 to aid in the systematic intro- 
duction of Si-Ro-Set to manufacturers 
of woolen and worsted garments in 
the United States. 

Si-Ro-Set involves spraying gar- 
ments with a certified Si-Ro-Set solu- 
tion which is similar to one used in the 
“home permanent” wave. It imparts 
durable pleats or creases to wool ap- 
parel without altering the hand or life 
of the fabrics. Garments treated with 
the process will hold the crease through 
rain, high humidity, dry cleaning or 
total immersion. 

Mr. Schmitt said that chemical 
manufacturers signing agreements 
with the Commonwealth Scientific and 
Industrial Research Organization for 
use of the certification mark Si-Ro-Set 
may do so without charge. Garment 
manufacturers purchasing Si-Ro-Set 
concentrate will be automatically 
given a label license to use the concen- 
trate in accordance with the process or 
processes covered by any patents or 
patent applications owned by the 
CSIRO, without charge. 


Dr. Laxer 


Unique Shelter for Use with Jupiter 
Missile.—A shelter that opens like the 
petals of a flower has been devel- 


oped by the U. S. Army Engineer 
Research and Development Labora- 
tories, Fort Belvoir, Virginia, for use 
with the Jupiter intermediate-range 
ballistic missile. The weatherproof 
shelter is designed to cover the lower 
portion of the rocket, creating a water- 
tight seal that protects instruments 
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and equipment until the rocket is 
ready to be fired. 

The shelter has twelve 30-ft. alumi- 
num panels which are shaped like 
the petals of a flower. After installa- 
tion, they radiate from a circle around 
the missile. When preliminary ad- 
justments on the missile are com- 
pleted, synchronized motors raise the 
panels, closing them around the rocket. 
When the missile is ready for firing, 
electronic equipment takes over, causing 
the panels to open up and return to 
the ground. 

Doors in the 30-ft. sections allow 
workers to go inside for inspection 
of the launching apparatus. Other 
openings carry fuel lines to the missile. 
The interior of the shelter is coated 
with a paint capable of resisting tem- 
peratures of over 1000 F. 

The engineering firm of Barnes and 
Reinecke, Incorporated, Chicago, IIli- 
nois, holds the principal Army con- 
tract for production of the shelter. 


Auto-Optic Recording Pyrometer. 
An optical recording pyrometer that 
can measure extremely high tempera- 
tures with milliseconds response time 
has been announced by Avco Research 
and Advanced Development Division 
of Wilmington, Mass. Developed to 
measure the high surface temperatures 
of sample materials exposed to Avco’s 
10,000 C. electric-arc air-heater, the 
new pyrometer may be operated un- 
attended when conditions are hazard- 
ous to personnel. 

Even under adverse conditions such 
as an electric-arc experiment, the py- 
rometer has given reliable measure- 
ments of sample brightness tempera- 
tures within +20 C. It accurately 
measures surface changes in tempera- 
ture over an area of 3-in. diameter in 
the range of 1400 C. to 3000 C. A 
recorder may be plugged into the py- 
rometer to provide a complete record 
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of temperature versus time for the 
entire operation. 

Called the Auto-Optic Recording 
Pyrometer, the radiation-sensitive ele- 
ment is a high-vacuum phototube. 
Monochromatic light of wavelength 
0.65 micron is focused upon an opaque 
plate having a small opening which 
admits the radiant energy to the pho- 
totube. The phototube output is am- 
plified by a highly stable amplifier 
mounted in the pyrometer housing. 
The output may be read on the built- 
in meter or from an external recorder. 

A remote recorder attached to the 
pyrometer is used by Avco to study 
high-temperature properties of new 
materials being developed for possible 
use on space vehicles. 

The instrument is calibrated by 
comparing the output with the known 
brightness temperature of a standard 
tungsten lamp. 


TASI Will Double the Capacity of 
Submarine Telephone Cables.—A new 
system which will increase the number 
of telephone conversations carried by 
underseas telephone cables perhaps by 
as much as two times is under develop- 
ment at Bell Telephone Laboratories. 
The system is called TASI, which 
stands for Time Assignment Speech 
Interpolation. 

In a normal two-way telephone con- 
versation, one person on the average 
will be speaking only half the time, 
and even while he is speaking, there 
will be significant gaps and pauses in 
his speech. Thus, if the two direc- 
tions of transmission are separated, 
each transmission path is idle, on the 
average, more than half the time. If 
two conversations are interlaced to 
take advantage of these gaps, greater 
use can be made of existing transmission 
facilities. It is this interlacing process 
which leads to the designation Time 
Assignment Speech Interpolation. 
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In practice, the system would not 
work with only two talkers on one 
line since they would frequently be 
speaking at the same time. However, 
where a larger number of channels is 
available, such as a submarine tele- 
phone cable, an averaging effect oc- 
curs so that at any instant there is a 
greater probability of sufficient “free 
time” being available to accommo- 
date the larger number of conversa- 
tions. Increasing the capacity of a 
36-channel system to 72 channels is 
more feasible with TASI than dou- 
bling the capacity of a 5-channel sys- 
tem. 
TASI is essentially a group of high- 
speed switches. If, for example, 36 
cable channels were available, 72 
talkers could be connected. When 
there are more talkers than channels, 
the equipment will connect talkers who 
become active by disconnecting talkers 
who are silent at that moment. In 
turn, this disconnected talker will be 
assigned another momentarily inactive 
channel when he starts to speak again. 
A talker will be disconnected only 
when he is silent. 

When a talker starts talking, his 
voice actuates a speech detector. The 
speech detectors are scanned by a 
control circuit similar to a modern 
digital computer. When a talker be- 
comes “‘active’’ this control circuit 
initiates a coded tone burst consisting 
of a group of four audio tones which 
precedes the voice over an available 
cable channel. After the tone burst, 
which lasts only 10 milliseconds, the 
control circuit connects the talker to 
the same channel. The coded tones 
operate switches to connect the talker 
to the proper line at the receiving end. 
The tones are not heard since the 
listener is not connected while they 
are being transmitted. Whena talker 
is not ‘‘active’’ and his channel is 
needed, another coded tone burst is 
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transmitted over a separate signaling 
channel and severs the connection. 

Switching of talkspurts from one 
channel to another is accomplished in 
a few milliseconds by a time-division 
switch. Speech is sampled for about 
two microseconds and the resulting 
pulses steered to the appropriate idle 
channel during the sampling by the se- 
lective operation of transistor “‘gates’’ 
in each channel. This short interval 
of sampling makes it possible to sam- 
ple all active talkers 8000 times a 
second and then reconstruct the speech 
from the samples before it is trans- 
mitted over the undersea cable. 

The signaling system in TASI keeps 
the receiving end informed of the con- 
nections that the transmitting end has 
established. Four groups of audio 
tones are employed for signaling pur- 
poses—four tones in each of three 
groups and three in the fourth group. 
Each signal comprises one and only 
one tone from each group—if more or 
less are present, an error is indicated 
and appropriate steps taken to cor- 
rect it. ‘‘Connect”’ signals which pre- 
cede the voice signal at the beginning 
of the speaker's talkspurt are sent 
over the same channel as the speech. 
Disconnect and connection checking 
signals are sent over a separate chan- 
nel used only for this purpose. 

All of the circuits for TASI are 
completely transistorized. Terminals 
for doubling the capacity of the Trans- 
atlantic cable will require several thou- 
sand transistors of four different types, 
and tens of thousands of semiconduc- 
tor diodes and passive components. 


Microscope Color TV System. 
Specialized medical television equip- 
ment, including a highly sensitive 
camera tube that can see in the dark 
and a compact transistorized color TV 
pickup system developed by the Radio 
Corporation of America, was demon- 
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strated recently to leaders in medical 
education, at a two-day conference on 
television uses in post-graduate medi- 
caleducation. The sessions were held 
at the Clinical Center of the National 
Institutes of Health, under joint spon- 
sorship of the American Academy of 
General Practice and the Institute for 
the Advancement of Medical Com- 
munication, Inc. 

The demonstration program was ar- 
ranged by Dr. Murray Brown, Chief 
of Clinical Professional Education at 
the National Institutes of Health, who 
has directed experiments at the In- 
stitutes with various television tech- 
niques. 

Demonstrations included a showing 
of results obtained by scientists with 
television pickup at very low light 
levels by color and black-and-white 
systems employing an RCA-developed 
wide-space Image Orthicon camera 
tube. The tube is so sensitive it can 
view scenes in surroundings which 
appear entirely dark to the human 
eye. 


One example shown was a color film 
recording of the retina of the human 
eye, achieved with a system employing 
the wide-space Image Orthicon to per- 
mit a light level low enough to avoid 


disturbing the patient. In another 
application, the extremely sensitive 
tube was employed to view a demon- 
stration of intricate apparatus, illumi- 
nated only by an ordinary 75-watt 
light bulb, showing application of such 
pickup equipment in normal surround- 
ings without the high-level lighting 
ordinarily required for studio-type 
television pickup systems. The sub- 
ject was an animal specimen within a 
germ-free tank, and the televised im- 
ages were displayed by a projection 
system onto a large screen adjacent 
to the tank, making them clearly visi- 
ble to the audience. 

The conferees also were shown a 
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fully transistorized miniature color 
television pickup system developed at 
the RCA Laboratories by a research 
group under Dr. V. K. Zworykin, 
Honorary Vice President of RCA. 
The 20-lb. camera, employing devel- 
opmental 3-in. Vidicon pickup tubes, 
picked up color images of subjects un- 
der a light microscope for display on 
standard color television receivers. 
The experimental system, contained 
in two compact units weighing only 
65 lb is designed for closed-circuit use 
in industry, education, defense and 
research. 

Dr. Zworykin pointed out that the 
system includes equipment enabling 
it to produce a standard color televi- 
sion broadcast signal permitting its 
possible use in picking up scenes for 
transmission over the air, as_ well 
as over closed-circuit systems. He 
added that the equipment was to be 
used in April by Smith, Kline and 
French Laboratories, Philadelphia, to 
transmit color microscope images to 
five medical schools in the Philadelphia 
area during a closed-circuit television 
lecture series. 


Stereo Sound Broadcasting Now 
Compatible with Single Channel Re- 
ception.—Several commercial radio 
and television stations have recently 
begun offering stereophonic sound pro- 
grams, by broadcasting over two sepa- 
rate sound channels. In various ex- 
perimental arrangements, the two 
channels required have been selected 
from different combinations of AM, 
FM, and TV channels. The listener 
then spaces the appropriate receivers 
properly in his home. Results have 
been sufficiently favorable that more 
and more broadcasters are now con- 
sidering offering stereo programming. 

The major obstacle to vastly in- 
creased use of this type of stereo 
broadcasting, however, is the person 
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who listens with only one receiver. 
If the broadcaster tries for the full 
stereo effect, the sound the single 
channel listener hears comes from only 
one of the two widely-spaced micro- 
phone pickups, and he misses a por- 
tion of the program. What he does 
receive is poorly balanced, because 
of the placement of the pickups in re- 
lation to the sound sources. The 
broadcaster has had to dilute the 
stereo effect in order to preserve satis- 
factory reception for the single chan- 
nel listener. 

Now this single channel problem 
may be eliminated without affecting 
the stereo listener, through the use of 
a new “compatibility” circuit which 
has been developed by F. K. Becker 
of Bell Telephone Laboratories. The 
circuit depends for success on a psy- 
cho-acoustic phenomenon known as 
the “Precedence Effect.’’ This ef- 
fect operates in such a manner that 
when a single sound is reproduced 
through two separate loudspeakers, 
but is delayed several milliseconds 


(thousandths of a second) in one, the 
listener will ‘“hear’’ the sound as if 
it came only from the speaker from 


which he heard it first. He will judge 
the second loudspeaker to be silent. 
The amount of the Precedence Effect 
depends somewhat on the length of 
delay built-in; at 10 milliseconds, the 
sense of direction for an average ob- 
server is as if the echo were 8 to 10 
decibels softer than the sound preced- 
ing it. 

In the new development, the cir- 
cuits between the microphone pick- 
ups and their corresponding radio or 
TV transmitters are cross connected 
through two delay lines, each with 
its own buffer amplifier. Because of 
these cross connections, music or voice 
signals from the left microphone are 
transmitted directly to the left loud- 
speaker in the listener’s home, while 
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the same signal is slightly delayed 
before reaching the speaker to his 
right. The stereo listener will hear 
the sound as if it came only from the 
left loudspeaker because of the Preced- 
ence Effect. Conversely, the sound 
from the right microphone goes direct 
to the right speaker, but is delayed 
before reaching the left speaker, and 
is therefore unheard. Thus, the brain 
of the stereo listener localizes the 
sound he hears as coming direct from 
each of his two speakers, and full 
stereophonic effect is maintained. 

However, monophonic reception is 
completely compatible with this, since 
a listener to each single channel hears 
the total sound from both micro- 
phones in a balanced reproduction. 
The slight delay built in apparently 
does not affect his reception at all, 
according to subjective tests performed 
at Bell Laboratories. 

Typical operating conditions in lis- 
tening tests conducted today had a 
time delay of 10 milliseconds, with the 
volume of the delayed channel 13 db 
softer than that of the direct. (The 
Precedence Effect, which was first 
discovered in 1933, operates over a 
delay range of from 5 to 35 milli- 
seconds. ) 

This development should make it 
possible for many more broadcasters 
to offer double-channel stereo pro- 
gramming without diluting the stereo 
effect or penalizing the single channel 
listener, who will make up the ma- 
jority of their audience. 


Solar Cookers for Underprivileged. 
University of Wisconsin anthropolo- 
gists are stationed in the Arizona desert 
and in central Mexico to watch In- 
dians bake their tortillas on solar 
cookers developed by the engineers 
in the UW Solar Energy Laboratory. 
They hope to find out for the engi- 
neers if the cooker, carefully built and 
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tested in the laboratory, is effective in 
constant daily use; and if it is, whether 
the people for whom it was developed, 
who live in areas where only sunlight 
and poverty are plentiful, will be 
willing to change their customs and 
use it every day. 

Four Indian tribes in Arizona have 
been cooperating with the anthropolo- 
gists in the tests since last summer, 
and during the second semester two 
selected villages in Mexico will also 
be under constant observation. 

Although scores of the cookers have 
been distributed in various under- 
developed areas of the world, no at- 
tempts have been made until now at 
extensive day-by-day studies of their 
efficiency and social acceptability, ac- 
cording to Prof. Milton Barnett, who 
heads the project for the anthropolo- 
gists. 

If the tests are successful, both 
scholars and Indians will profit: the 
anthropologists in valuable material 
on the hows and whys of cultural 
change, the engineers in the satisfac- 
tion of building something of benefit to 
under-privileged mankind, and the In- 
dians—Mohaves, Chemehuevi, Hopi, 
and Navajo in Arizona, and the Otomi 
in Mexico—with a cheap, quick 
method of cooking their food. 

The anthropologists spend a prelim- 
inary period with the Indians to gather 
basic ethnographic data on their phys- 
ical characteristics and social and cul- 
tural traits and to establish friendly 
relations with the families who will 
participate in the testing. Only then 
do they bring in the saucer-shaped and 
plastic-coated cookers which collect 
and concentrate the rays of the sun on 
a pot or pan mounted on an attached 
rack. 

The heat is intense enough to boil 
a quart of water in 10 minutes. They 
carefully demonstrate the cookers un- 
til the natives have mastered the 
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operation, and then they observe per- 
sonal and cultural reactions to them. 

“The men in Engineering are-on 
their toes,’ Barnett says. ‘Because 
one of the Hopi women suggested that 
hot water for laundry was a problem, 
they designed a solar unit to heat 
water. They are working on a design 
for a do-it-yourself kit to cut down 
the cost of the heater and make it 
more available to people short of cash ; 
and two students in chemical engi- 
neering are designing a field sterilizer 
for doctors and nurses to use in iso- 
lated places. 

The UW Solar Energy Laboratory 
is one of the few in the world carrying 
on a broad research program in use of 
sunlight as a source of energy. The 
program was started in 1955 under 
direction of Prof. Farrington Daniels, 
chairman of the chemistry department 
and internationally recognized author- 
ity on solar energy. 


New Error-Correction Code for 
Bursts of Errors.—Lightning flashes 
and other electrical disturbances which 
cause static and noise on communica- 
tion lines may result in groups or 
bursts of errors in the data being 
transmitted over these lines. These 
errors can be largely eliminated by a 
new error-correcting code developed 
by Dr. D. W. Hagelbarger of Bell 
Telephone Laboratories. The termi- 
nal equipment required for this new 
code is simple and inexpensive, and 
synchronization is relatively easy to 
maintain. 

Previous error detecting and cor- 
recting codes either could not handle 
adjacent errors, or required compli- 
cated terminal equipment and _ pre- 
sented difficult synchronization prob- 
lems. The new code is applicable to 
systems where the data must be ac- 
cepted and delivered continuously, 
rather than in batches, according to 
the inventor. 
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Error correcting systems are not 
required in the transmission of speech 
or television signals, because the in- 
formation is basically analog in char- 
acter, and high degrees of redundancy 
exist. This reduces the detrimental 
effect of any errors which may occur. 
However, in telegraph lines, or in data 
transmission such as the Dataphone 
system recently announced by Bell 
Laboratories, information trans- 
mitted in a discrete pulse-no pulse 
form, and mutilation of a single digit 
could throw a complete set of data in 
doubt. In this type of transmission, 
bursts of errors could be critical, and 
correction methods become an essen- 
tial part of transmission equipment. 

Terminal equipment using the new 
code can be designed to handle prac- 
tically any length error burst which 
system analysis indicates is required. 
In general, the shorter the maximum 
burst length to be corrected, the smaller 
and simpler will be the terminal equip- 
ment. Also, a short burst length will 
result in less guard space, or “clean 
data section,’ which must follow the 
burst before another group of errors 
can be corrected. 

In its simplest form, the coding 
system uses alternate data digits and 
check digits, giving a redundancy of 
one-half. If such a system is designed 
to correct error bursts of length six 
or less (three data digits and three 
check digits), the encoder consists 
basically of a shift register of length 
seven. The data digits enter the first 
position, and are shifted through the 
register before being transmitted. At 
each shift, a check digit is computed 
which makes parity (the sum of 1’s) 
of the check digit and the data digits 
in the 1st and 4th positions even (zero 
or two). This check digit is trans- 
mitted soon after it is computed, pre- 
ceding the transmission of its nearest 
associated message digit by seven 
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digits. Data and check digits emerge 
from the coding system alternately, 
forming the coded message. 

At the decoder, the received mes- 
sage is separated into check and data 
digits, each group entering a separate 
shift register. There are two copies 
of the parity circuits: one checks the 
parity relationship among data digit 
1, data digit 4, and checking digit 7, 
while the second checks the parity 
among data digits 4 and 7 and check 
digit 10. The decoding rule is: when- 
ever both parity circuits fail (by indi- 
cating an uneven parity), the data 
digit in position 4 is changed while 
being shifted to position 5. Ifonly one 
parity check fails, no change is made. 

With the redundancy of 4, the syn- 
chronization problem is minimized, 
since the data digits and check digits 
alternate. 

After any burst of length 6 or less, 
a 19-digit errorless message or guard 
space is enough to fill the decoder 
shift registers completely and prepare 
the decoder for another burst. 

The redundancy of these codes can 
be reduced if desired, or they can be 
used to correct for bursts longer than 
6, without difficulty, but with simul- 
taneous increases in terminal equip- 
ment costs and in the guard space 
required. 

The decoder can be equipped to 
detect bursts longer than those it is 
designed to correct. For instance, a 
code to correct error bursts of length 
six will detect all but seven of the 512 
possible bursts of length 9 or less, and 
operate a warning device. 

Evaluation of the new code is now 
being conducted under transmission 
line conditions by Dr. F. E. Froehlich 
of Bell Laboratories. Initial results 
indicate substantial improvement fac- 
tors in code redundancies of 4 and }. 
The evaluation program is continuing. 
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The International Yard and Pound. 

The Directors of the foliowing stand- 
ards laboratories—Applied Physics Di- 
vision, National Research Council, 
Ottawa (Canada), Dominion Physical 
Laboratory, Lower Hutt (New Zea- 
land), National Bureau of Standards, 
Washington (United States of Amer- 
ica), National Physical Laboratory, 
Teddington (United Kingdom), Na- 
tional Physical Research Laboratory, 
Pretoria (South Africa), National 
Standards Laboratory, Sydney (Aus- 
tralia)—-have discussed the existing 
differences between the values assigned 
to the yard and to the pound in differ- 
ent countries. To secure identical 
values for each of these units in precise 
measurements for science and _ tech- 
nology, it has been agreed to adopt an 
international yard and an international 
pound having the following definition : 
the international yard equals 0.9144 
meter; the international pound equals 
0.453 592 37 kilogram. 

It has also been agreed that, unless 
otherwise required, all nonmetric cali- 
brations carried out by the above labo- 
ratories for science and technology on 
and after July 1, 1959, will be made in 
terms of the international units as de- 
fined above or their multiples or sub- 
multiples. 

It will be noted that the Interna- 
tional Inch is approximately 2 parts 
per million shorter than the inch pres- 
ently used by the National Bureau of 
Standards, and somewhat less than 2 
parts per million longer than the inch 
now used by the National Physical 
Laboratory. ‘To avoid possible con- 
fusion, during the transition period 
National Bureau of Standards calibra- 
tions of length or mass expressed in 
English units will embody a statement 
indicating clearly the unit which has 
been used if the choice introduces a 
significant difference in the calibration 
values. Furthermore, if the accuracy 
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of the calibration is such that the certi- 
fied values would be the same in either 
“International” units or the older 
units, the qualifying adjective “Inter- 
national” will not be used, that is, the 
values will be expressed, for example, 
as so many inches or pounds. 

This unit shall be referred to as the 
American Survey Foot. Inasmuch as 
there is little or no interchange of sur- 
vey data, where the foot measurements 
are used, with industrial and scientific 
data, where the international units will 
be used, it is anticipated that no con- 
fusion will result from this dual usage. 
For example, base line surveys which 
might enter into a velocity of light de- 
termination would invariably be made 
in terms of meters. 

The values of the pounds currently 
in use in the United States, United 
Kingdom and Canada are as follows: 
1 U.S. pound = 0.453 592 4277 kg; 1 
British pound = 0.453 592 338 kg; 1 
Canadian pound = 0.453 592 43 kg; 1 
International pound = 0.453 592 37 kg. 

The relative differences in the vari- 
ous pounds are substantially less than 
in the yards but since masses can be 
measured with greater accuracy than 
lengths, the differences can be signifi- 
cant. The present British pound is 
about one part in ten million smaller 
than the international pound, whereas 
the U. S. and Canadian pounds are 
about one and one-half parts in ten 
million larger. 

The conversion factor for the inter- 
national pound was selected so as to be 
exactly divisible by 7 to give the fol- 
lowing value for the grain: 1 Interna- 
tional grain = 0.06479891 gram. The 
grain is the common unit in avoirdu- 
pois, apothecary and troy pounds. 
There are 7,000 grains in the avoirdu- 
pois pound, and 5,760 grains in both 
the apothecary pound and the troy 
pound. “ 
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